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I .  INTRODUCTION 


This  thesis  deals  with  satellite  tracking  of  relatively 
low  altitude  targets.  The  long  term  goal  of  this  project  is 
to  determine  an  effective  general-coverage  satellite  orbital 
pattern.  The  short  term  goal  of  this  thesis  is  to  explore 
the  possible  use  of  various  coordinate  systems  and  to  design 
a  satellite  observer.  A  target  such  as  a  slow  moving 
aircraft  is  assumed  for  the  basic  development. 

Chapter  2  contains  an  overview  of  basic  satellite 
mechanics  and  dynamics  as  well  as  a  brief  discussion  of 
satellite  detection  and  tracking  equipment.  Equations 
regarding  orbital  motion  will  also  be  introduced.  Also 
chapter  2  provides  an  overview  of  the  entire  thesis. 

The  theory  and  mathematics  of  satellite  and  target 
motion  for  a  specific  case  are  developed  in  Chapter  3.  The 
geometry  of  motion  on  a  spherical  surface  is  detailed  in 
that  chapter.  Specifics  such  as  geometric  swath  width. 

1 ine-of -sight  and  precession  are  furnished.  This  leads  into 
an  in-deapth  examination  of  coordinate  systems.  Various 
attempts  to  find  a  suitable  set  of  reference  axes  for 
analytic  studies  of  the  tracking  problem  are  presented.  It 
is  desired  to  achieve  a  target  and  observation  model  that 
has  a  convenient  basis  for  analysis  such  as  either  linear  or 


bilinear  in  state  space.  Several  trial  coordinate  systems 
were  developed  mathematically  but  found  unsuitable.  One 
trial  system  is  considered  satisfactory  under  certain 
assumptions  and  is,  therefore,  completely  developed. 

Chapter  4  presents  a  discourse  on  observability  in 
general,  and  non-linear  observability  in  particular.  The 
concept  of  connectedness  and  univalence  is  discussed.  The 
chosen  system  is  then  analyzed  for  observability. 

The  topic  of  observers  is  presented  in  Chapter  5.  A 
background  on  basic  obsevers  is  provided.  Finally  an 
observer  is  designed  and  simulated  by  computer  for  the 
prefered  system. 

Chapter  6  offers  conclusions  and  recommendations  for 


further  studies 


THE  SATELLITE  IN  ORBIT 


Before  an  understanding  of  satellite  tracking  can  be 
achieved  it  is  essential  to  have  an  understanding  of  basic 
satellite  dynamics  and  satellite  sensors.  The  concepts  and 
equations  that  govern  the  motion  of  a  satellite  in  orbit 
around  the  earth  are  based  on  the  physical  laws  put  forth  by 
kepler  and  Newton.  These  laws  of  motion  apply  to  artificial 
satellites  as  well  as  planets  and  moons.  The  physical 
geometries  and  forces  are  the  same.  The  following 
explanations  and  derivations  form  a  basis  for  orbital  motion 
of  a  satellite. 

A.  THE  SHAPE  OF  THE  EARTH 

In  reality  the  earth  is  neither  exactly  round  nor 
exactly  a  rigid  body.  The  constituents  can  be  molded  and 
deformed  to  a  certain  degree.  The  earth's  rotation  causes 
it's  own  materials  to  feel  a  force  that  pulls  tangentially 
into  space.  Fortunately  the  earth's  materials  also  undergo 
a  constant  centripetal  (center  seeking)  acceleration  to  keep 
them  in  place.  This  centripetal  acceleration,  which  is 
directed  towards  the  center  of  curvature,  is  provided  by 
earth's  gravity. 

The  gravitational  force  produced  by  the  earth  is 
directed  towards  the  center  of  the  earth.  However,  the 
momentum  produced  by  the  earth's  rotation  pulls  at  the 


earth'3  materials.  Each  particle  of  earth  has  both  these 
forces  acting  on  it.  The  earth's  surface  at  the  equator  is 
actually  moving  faster  around  the  rotation  axis  than  the 
surface  material  at  or  nearer  the  poles.  Figure  2.1 
illustrates  this.  Therefore  the  materials  on  the  surface  at 
the  equator  feel  this  tangential  pull  the  most.  Analysis 
shows  that  the  result  is  an  earth  shaped  like  an  oblate 
ellipsoid . 

An  oblate  ellipsoid  can  be  described  as  a  sphere  which 
has  been  compressed  along  the  polar  axis  and  therefore 
bulges  at  the  equator.  A  cross  sectional  slice  through  the 
poles  yields  an  ellipse.  A  slice  along  the  equatorial  plane 
yields  a  circle. 

This  oblateness  (measure  of  the  earth's  flattening)  of 
the  earth  is  only  slight.  The  equatorial  diameter  is 
calculated  to  be  12,757  km  whereas  the  polar  diameter  is 
calculated  to  be  12,714  km.  This  is  only  a  43  km  difference 
which  is  about-  1  part  in  297. 

B.  GREAT  AND  SMALL  CIRCLES 

The  earth'3  equatorial  bulge  can  be  temporarily 
disregarded  in  order  to  consider  the  earth  as  a  sphere 
turning  beneath  an  orbiting  satellite.  The  intersection  of 
a  plane  passed  through  the  center  of  the  sphere  and  the 
sphere  is  the  largest  circle  that  can  be  drawn  on  the 
sphere's  surface.  This  is  known  as  a  great  circle.  The 
shortest  distance  between  any  two  points  on  the  surface  of 


the  earth  is  an  arc  of  a  great  circle.  The  earth's  equator 
is  an  example  of  a  great  circle.  There  are  an  infinite 
number  of  approximate  great  circles  on  the  earth's  surface 
since  there  are  an  infinite  number  of  planes  that  can  be 
passed  through  the  earth's  center. 

There  is  a  series  of  24  special  great  circles  that  pass 
through  the  earth's  poles  and  are  evenly  spaced  from  each 
other.  These  are  called  meridians  and  they  intersect  with 
the  equator  at  right  angles.  Meridians  join  at  both  poles. 
Forming  other  right  angles  to  the  Meridians  are  parallels. 
Parallels  are  small  circles  created  by  passing  planes 
through  the  earth  parallel  to  the  equatorial  great  circle. 
See  Figure  2.2. 

In  order  to  determine  precise  locations  on  the  earth's 
surface,  latitudes  and  longitudes  are  utilized.  The 
longitude  of  a  specific  place  refers  to  the  arc  (measured  in 
degrees)  of  a  parallel  between  that  place  and  the  prime 
meridian  (which  passes  through  Greenwich,  England  for 
reasons  of  history).  Longitudes  run  east  and  west. 
Latitudes,  however,  run  north  and  south.  The  latitude  of  a 
specific  place  may  be  defined  as  the  arc  (in  degrees)  of  a 
meridian  between  that  place  and  the  equator.  Figure  2.3 
shows  an  example.  When  the  earth  is  considered  as  having  an 
ellipse  as  a  cross  section  (which  it  has)  instead  of  a 
circular  one,  the  length  of  a  degree  of  latitude  is  slightly 
greater  at  the  poles  than  at  the  equator. 


F  i  gure  ■£'.  3 
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C.  ORBITAL  MOTION 

A  presentation  concerning  orbital  motion  is  basic  to 
satellite  tracking.  Gravitational  attraction  and  momentum 
combine  to  keep  a  satellite  in  perpetual  orbit  above  the 
earth  (in  the  absence  of  friction  and  external  celestial 
attractions) . 

1 .  Gravity 

The  understanding  of  gravity  and  orbital  motion  can 
be  traced  back  to  Johannes  Kepler  (1571-1630)  and  Isaac 
Newton  (1642-1727) .  Kepler  is  known  predominantly  for  his 
three  laws  of  planetary  motion.  The  laws  are  as  follows; 

(1)  Each  planet  moves  about  the  sun  in  an  orbit  that  is  an 
ellipse,  with  the  sun  at  one  of  the  foci  of  the  ellipse. 

(2)  The  straight  line  joining  a  planet  and  the  sun  sweeps 
out  equal  areas  in  space  in  equal  intervals  of  time. 

(3)  The  squares  of  the  sidereal  periods  of  the  planets  are 
in  direct  proportion  to  the  cubes  of  the  semimajor  axes  of 
their  orbits. 

The  above  three  laws  apply  equally  well  for 
artificial  satellites  as  they  do  for  planets.  Newton 
restated  and  clarified  Kepler's  laws.  Newton  was  an 
advocate  of  rigorous  proofs  whereas  Kepler  preferred  to 
state  empirical  laws  based  on  observations. 


Kepler's  first  law  stated  mathematically  is 


where:  f “centripetal  force  needed  for  circular  orbit, 

*=*as*  of  planet  (or  satellite)  in  orbit, 
v  =  velocity  of  orbiting  body, 

r  3  diatance  between  sun  and  planet  (or  satellite  and 
earth) , 

Kepler's  third  law  restated  mathematically  is 

(m,  +  nOp*  s  4xJeVG  (2.2) 

where  ml  and  m2  refer  to  the  masses  of  two  bodies  that 
revolve  mutually  about  each  other, 
p  3  period  of  revolution, 
a  3  semimajor  axis  of  relative  orbit, 

G  *  universal  gravity  constant  (  6.67  x  10  ^ Vk<j  •  j “  ) 

Newton  also  explained  the  reasons  behind  kepler's 
original  observations.  That  is,  Newton  supplied  the  laws  of 
motion  which  are  at  the  root  of  classical  mechanics. 

These  laws  are  as  follows; 

(1)  Every  object  remains  at  rest  or  in  uniform  motion 
unless  an  external  force  acts  upon  it. 

(2)  The  product  of  the  mass  of  an  object  and  its 
acceleration  vary  directly  as  the  resultant  force,  and 
the  change  in  motion  takes  place  in  the  direction  of 
that  force. 

(3)  For  every  action  there  is  an  equal  (in  magnitude)  and 
opposite  (in  direction)  reaction. 

In  the  case  of  orbiting  planets  or  satellites, 

gravity  is  the  unseen  force  referred  to  in  Newton's  first 

law.  Newton  also  postulated  that  every  particle  of  matter 

in  the  universe  attracts  every  other  particle.  The  force  of 


this  attraction  is  proportional  to  the  product  of  their 
masses  and  inversely  proportional  to  the  square  of  the 
distance  which  separates  them. 

In  equation  form  this  is 

F  =  Gm,  Ml  / d 1  (2.3) 

where  d  is  the  distance  between  the  center  of  mass  1  and 
mass  2. 

Orbital  motion  can  be  understood  in  terms  of 
Newton's  laws  previously  stated.  The  critical  factors  in 
putting  a  satellite  in  orbit  (by  the  traditional  launch 
method)  are  speed  and  direction  of  movement  at  burnout. 
Burnout  is  when  the  rocket  engine  shuts  off  and  the 
satellite  behaves  as  an  astronomical  object. 

As  a  satellite  follows  it's  orbital  path  it 
continuously  falls  toward  the  earth  due  to  the  earth's 
gravitational  pull.  However,  the  satellite's  momentum 
prevents  it  from  really  being  pulled  any  nearer  to  the 
earth.  The  satellite's  orbit  is  the  result  of  two  main 
forces.  The  momentum  of  the  satellite  is  a  measure  of  its 
state  of  motion.  The  inertia  of  the  satellite  (recall 
Newtons  first  law)  is  that  property  that  causes  the 
satellite  to  resist  acceleration  and  travel  in  a  straight 
line.  For  the  satellite  to  move  in  a  circular  path  rather 
than  in  a  straight  line,  it  must  continually  suffer  an 
acceleration  toward  the  center  of  the  circle.  This 
acceleration  is  centripetal  acceleration.  The  central  force 


that  produces  the  centripetal  acceleration  is  that 
gravitational  attraction  between  the  satellite  and  the 
earth.  Figure  2.4  illustrates  this  in  vector  form. 

The  period  of  an  artificial  satellite  is  directly 
related  to  the  size  of  it's  orbit;  i.e,  the  closer  the 
satellite  is  to  the  earth  the  faster  it  travels. 

2 .  The  Center  of  Hass 

The  center  of  mass  can  be  defined  as  that  point 
within  a  system  that  either  remains  fixed  or  moves  as  if  the 
entire  mass  of  the  system  were  concentrated  at  that  point. 
For  a  number  of  particles  the  center  of  mass  is  defined  as 

£cn  *  >.?£?.  (2.4) 

An  example  of  this  for  a  two  body  system  is 

Xc*v  *  *jjc,  (2.5) 

m,  ♦ 

Hass  1  is  a  distance  x(l>  from  an  arbitrary  origin, 
and  mass  2  is  a  distance  x(2>  from  the  origin.  X(cm) 
represents  the  location  of  the  center  of  mass.  A  solid 
object  can  be  thought  of  as  a  collection  of  a  great  many 
particles.  In  this  case  equation  2.4  applies. 

Very  often  it  is  useful  in  certain  calculations  to 
treat  a  particular  body  as  if  its  entire  mass  is 
concentrated  at  one  point.  That  one  point  is  the  center  of 


mass 


It  can  be  shown  tRef.  1]  that  a  large  sphere 
attracts  other  bodies  as  though  the  sphere's  mass  is  all 
located  at  the  center.  This  holds  true  as  long  as  the 
sphere  is  of  uniform  density  or  is  made  up  of  concentric 
shells  each  being  of  uniform  density.  The  earth  is  often 
approximated  as  being  a  sphere  made  up  of  uniformly  dense 
shells  CRef .  23 . 

3 .  Equations  of  Motion 

A  satellite  in  orbit  around  the  earth  will  feel  by 
far  the  strongest  gravitational  pull  from  the  earth  itself. 

*i 

Other  astronomical  bodies,  such  as  the  moon  or  sun,  do  exert 
force  disturbances  but  are  too  far  away  to  exert  a  very 
strong  pull  on  the  satellite.  Assuming  the  gravitational 
acceleration  of  the  satellite  is  due  only  to  the  earth's 
gravitational  attraction,  g  can  be  defined  as 

g  =  Gm  .  /r  *  (2.6) 

The  value  of  g  does  vary  with  respect  to  altitude 
and  latitude  (since  the  earth  is  not  really  an  exact 
sphere).  Tables  1  and  2  depict  the  variations  of  g.  The 
values  of  g  are  close  enough  to  allow  an  approximation  of 
the  earth's  shape  as  spherical  for  most  appl icat ions . 

An  orbiting  satellite  is  considered  to  be  in  the 
earth's  gravitational  field.  The  vector  symbol  g  refers  to 
this  field  and  is  defined  at 
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TABLE  1 


VARIATION  OF  g  WITH  LATITUDE  AT  SEALEVEL 


Latitude 

g  ( m /s'*  ) 

0* (equator ) 

9.78039 

10* 

9.78195 

20° 

9.78641 

30* 

9.79329 

40° 

9.80171 

50* 

9.81071 

60  * 

9.81918 

70* 

9.82608 

80° 

9.83059 

90* 

9.83217 

TABLE  2 

VARIATION  OF  g  WITH  ALTITUDE  AT  45*  LATITUDE 


Altitude(km ) 

g ( m/s 1 

0 

9.806 

1 

9.803 

4 

9.794 

8 

9.782 

16 

9.757 

32 

9.710 

100 

9.600 

500 

8.530 

1000 

7.410 

380000 

0.003 

g 


F/m 


(2.7) 


When  an  object  is  a  distance  h  above  the  earth's 
surface  it  has  potential  energy  defined  as 

U  =  - ( -mg ) h  =  mgh  (2.8) 

Gravity  is  a  conservative  force  pointing  to  the 
earth's  center  and  has  the  value  <-mg).  The  gravitational 
force  exerted  on  an  object  in  this  field  can  be  derived  from 
the  potential  energy  equation  as 

F  =  -du  =  d  f-GMml  =  -GMm  (2.9) 

dr  dr  l  r  j  r2 

A  more  general  form  for  potential  energy  invoking 
the  universal  constant  G  is 

U  =  -GMm/r  (2.10) 

Potential  energy  can  be  converted  into  kenetic 
energy  CRef.  31.  As  an  object  falls  to  earth  it  loses  its 
potential  energy  as  its  height  decreases  but  gains  kenetic 
energy  as  its  velocity  increases.  Kenetic  energy  is 
mathematically  defined  as 

K  =  GMm/2r  =  mv‘/2  (2.11) 

The  total  mechanical  energy  is 

E  =  K  ♦  U  (2.12) 

Orbital  motion  can  be  considered  as  motion  in  a 
plane.  For  one  full  revolution  around  the  earth  the 
satellite  remains  in  a  plane  as  it  traces  out  an  ellipse  or 
circle . 
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The  velocity,  period  and  altitude  of  a  satellite  are 
all  inter-related.  The  period  of  a  satellite  [Ref.  23  can 
be  defined  as 

X-  2ir(R  *  h)/v  ,  (2.13) 

where  v  is  the  velocity  of  the  satellite. 

The  centripetal  force,  is  often  defined  as 

F  =  m3v,/(R  «•  h)  (2.14) 

where  m  is  the  mass  of  the  satellite. 

It  is  obvious  to  the  casual  observer  that  the 
centripetal  force  is  also  in  part  a  function  of  the 
satellite's  altitude  and  velocity.  Polar  coordinates  are 
used  in  figure  2.5  to  illustrate  the  elliptical  orbit 
variables. 

4 .  Uniform  Circular  Motion 

Sometimes  a  satellite's  elliptical  orbit  is  that 
perfect  ellipse  the  circle.  If  the  satellite's  orbital 
velocity  remains  unchanged  in  a  circular  path  then  the 
satellite  is  moving  with  uniform  circular  motion. 

Since  the  orbit  of  a  two  body  system  (earth  and 
satellite)  is  often  in  a  plane  it  is  possible  to  use  polar 
coordinates  to  develop  basic  equations.  Figure  2.6  shews 
the  relation  between  polar  and  rectangular  coordinates. 
Velocity  in  polar  coordinates  for  circular  motion  can  be 
expressed  as 
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(2.15) 


where  0  refers  to  motion  in  the  theta  (0)  direction. 
Acceleration  is  found  by  taking  the  derivative  of  the 
velocity 

a  =  -r0Jr  ♦  r08  (2.16) 

Energy  can  be  expressed  as 

E  =  K  U  =  m(rJ0')/2  -  GMm/r  (2.17) 

and  angular  momentum  as 

L  =  mr *  0  (2.18) 

Uniform  circular  motion  is  easier  to  deal  with  than 
non-uniform  circular  motion.  There  i3  only  an  angular 
component  of  velocity  present  in  the  equation.  It  is  often 
useful  to  simplify  the  satellite  dynamics  in  order  to  more 
easily  manipulate  equations. 

5 .  Equations  of  Motion  in  Polar  Coordinates 

Consider  the  special  case  of  a  satellite  in  a 
circular  orbit.  This  can  be  treated  as  motion  in  a  plane. 
Polar  coordinates  are  especially  usefull  for  this  case.  The 
two  variables  that  define  motion  and  position  are  the  radial 
coordinate,  r,  and  the  angular  coordinate,  0,  conversion 
between  rectangular  coordinates  and  polar  coordinates  is  as 
follows . 

r  =  \)  X1  *  Y*  X  =  rcos  ( 0 ) 

0  =  arctan(Y/X)  Y  =  rsin(0) 

Figure  2.7  illustrates  this. 
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Polar  coordinates  are  often  used  any  time 


curvilinear  motion  occurs.  In  Figure  2.8  basic  curvilinear 
motion  is  shown.  The  differential  length  dr  has  components 
in  both  the  r  and  8  directions.  The  vector  r  can  be 
expressed  as  rn  where  n  is  a  unit  vector  in  the  radial 
direction . 

Velocity  can  be  defined  as; 

v  *  dCr)  =  d(rn,~>  =  d(r*  nr)  ♦  rdn,  (2.19) 

dt  “cTE  dt  dt 

Rectangular  coordinates  are  used  once  again  to  achieve  an 
expression  for  the  time  derivative  of  the  unit  vector  n  and 

!Le  ;  - 

n_r  =  cos(0)i  ♦  sin(0)j  n^  =  -sin(8)i  *  cob  ( 8 )  ] 

nr  =  -sin(8)i8  ♦  cos(8)j8  3  <-sin(9)i  *  cos<9)3)0  3  She 
=  -cos(8)i9  -  sinCQij©  -  -(cos(8)i  ♦  sin(0)3)0  3  -0n, 

Velocity  can  now  be  expressed  as; 

V  =  rn  *  r8n 

—  r  —  <9 

Following  the  same  rules  of  differentiation  acceleration  can 
be  derived. 

a_  =  V  =  rn  r  *  r8n  +  r8n  ♦r9n(.,  -r8Jn 
Simplifying , 

a  =  (r  -  r0')n  *  (2r9  ♦  r8)n  (2.20) 

—  -  r  — 

For  circular  motion  the  radial  component,  r,  remains 
constant.  Only  0  varies.  Therefore; 

V  =  r 8n  y  (2.21) 

a 


-r8'n  *  r9n 


(2.22) 


igure  £.7 


Coordinate  Conversion 


F 1 gure 


Vector 


The  acceleration  term  directed  in  the  negative  radial 
direction  is  the  center  seeking  acceleration. 

Figure  2.9  depicts  the  directions  for  the  force  vectors. 
The  sum  of  forces  is  calculated  as  follows; 

ZF  =  ma  =  m  (r-r0  *  )  n  (2.23) 

-  r  —  r  — r 

ZF^  =  ma  =  m<r0  ♦  2r0)n  (2.24) 

-e  -e  -e 

Using  differential  equations  and  state  form  motion  can  be 
represented  as  in  equations  2.27  and  2.28. 


dr/dt  =  Vr 

<2.25) 

d0/dt 

(2.26) 

dV 

/dt 

=  Vg  /r  -  k/r  +  U  -  6,- 

<2. 27) 

d- 

/dt 

=  -2Vri<^/r  *-(U  -6^/r 

(2.28) 

where  U,  =  Fr/m  and  U^,  =  Fe/n.  Both  variables  represent 
thrusts.  The  symbols  6r  and  6e  are  used  to  represent 
disturbances  which  are  suumed  to  be  zero  in  this  case. 

D.  SATELLITE  SENSORS 

Various  sensors  on  board  the  satellite  are  employed  to 
keep  the  satellite  itself  on  track,  and  to  detect  and  track 
a  possible  target.  The  satellite  guidance  system  uses 
sensors  for  measuring  certain  vehicle  dynamic  variables. 
Such  variables  include  satellite  acceleration,  velocity, 
position  and  angular  velocity.  Guidance  systems  generally 
determine  these  variables,  compare  the  received  information 
with  desired  parameters,  and  generate  correction  commands, 
various  types  of  sensors  include  inertial  7  optical  and 
radio  sensors.  The  satellite  may  make  use  of  gyroscopes  to 


establish  a  set  of  reference  axes  on  board.  Active  thrust 
control  is  used  to  align  the  body  axes  with  the  reference 
axes.  Once  the  satellite  has  established  the  craft's 
position  and  attitude  active  control  components  are 
activated.  These  components  respond  either  to  a  telecommand 
from  a  ground  station  or  to  the  satellite's  on-board 
computer  in  order  to  correct  the  satellites's  status.  In 
some  applications  where  perturbations  are  likely,  feedback 
control  is  added  to  the  active  control  of  a  satellite.  In 
some  cases  optical  sensors  are  used  in  conjunction  with 
inertial  equipment.  Optical  sensors  include  such  devices  as 
sun  finders,  sun  sensors,  planet  sensors,  celestial 
trackers,  and  horizon  scanners.  The  gyroscopes  and 
accelerometers  are  examples  of  inertial  sensors.  In  certain 
applications  the  inertial  sensors  provide  short-term 
stability  and  optical  sensors  provide  long-term  stability. 

The  gyroscope  is  an  inertial  guidance  system  that  is 
often  used  for  ships,  aircraft  and  spacecraft.  A  gyroscope 
is  a  device  that  possesses  a  high  rate  of  spin  about  an  axis 
of  symmetry  that  has  freedom  of  angular  rotation.  Examples 
of  gyroscopes  are  the  rate  gyro  and  the  integrating  gyro. 

Surveilance  satellites  normally  have  an  entirely 
different  set  of  sensors  for  target  acquisition.  Very  often 
this  type  of  satellite  carries  several  types  of  information 
gathering  devices.  High  resolution  optical  devices,  radar 
scanners  and  infrared  capacity  can  all  be  included  in  a 


satellite  payload.  Passive  systems  are  sensitive  receivers 
which  normally  provide  data  on  bearing  and  bearing  rate. 
These  systems  have  the  advantage  that  they  do  not  alert  the 
target  of  their  actions.  Optical  techniques  normally 
involve  lenses  and  cameras.  This  can  be  very  precise, 
especially  at  low  altitudes. 

Radar  is  an  example  of  an  active  sensor.  It  bounces 
radio  waves  off  a  chosen  target.  This  is  excellent  for 
determining  both  bearing  and  range.  Both  optical  and  radar 
methods  can  be  degraded  by  poor  atmospheric  conditions.  For 
the  radar  tracking  technique  noise  is  always  the  main 
limitation . 

Inverse  scattering  methods  are  being  developed  to  enable 
polarimetric  radar  to  obtain  better  target  information. 

This  procedure  involves  illuminating  the  target  with 
polarized  waves  and  observing  the  amplitudes  and  phases  of  a 
set  scattered  waves. 

Another  type  of  sensor  is  the  synthetic  aperture  radar 
system.  This  is  based  on  holography  methods.  Holography 
has  the  ability  to  record  three  dimensional  pictures  and 
focus  sharply  on  both  the  near  field  and  the  far  field  at 
the  same  time.  This  involves  recording  a  wave  interference 
pattern.  A  microwave  generator  is  used  to  provide  a 
constant  frequency  microwave  signal  as  well  as  a  reference 
wave.  This  is  considered  to  be  a  highly  precise  and 
accurate  sensor. 


Passive  sensors  are  normally  infrared  sensors.  There 
are  two  tajor  techniques  used  in  enfrared  detectors.  The 
two  main  types  are  called  photon  detectors  and  thermal 
detectors. 

In  photon  detectors,  the  technique  used  involves  a 
photon  of  infrared  radiation  being  absorbed  by  a  semi¬ 
conductor  electron  which  raises  its  energy  level  into  a 
conduction  band.  Photons  with  less  energy  than  the  band  gap 
produce  bo  signal  and  are  thereby  effectively  filtered  out. 

Thermal  detectors  function  be  sensing  the  temperature 
change  resulting  from  absorption  of  infrared  radiation  by  a 
suitable  element.  Usually,  this  absorbing  element  has  some 
temperature  sensitive  electrical  property  such  as 
resistivity  so  that  the  temperature  change  is  sensed 
electrically . 

One  device  often  used  in  passive  detectors  is  the 
radiometer.  This  device  is  a  broadband,  dual  frequency  , 
low  noise,  solid  state,  remote  controlled  mechansim. 
Radiometers  are  used  to  obtain  high  resolution  imagery  from 
low  altitude  satellite's. 

A  promising  type  of  detector  is  the  Silicide  Schottky 
diode  based  infrared  camera.  These  are  easy  to  manufacture 
and  have  excellent  performance.  The  camera  focal  plane 
consists  of  a  two  dimensional  array  of  metal  electrodes 
fabricated  on  a  silicon  substrate.  The  focal  plane  is  back 
illuminated.  When  an  infrared  picture  is  observed,  hot 


31 


carriers  are  emitted  from  the  focal  plane  electrodes.  An 
electronix  image  of  the  scene  is  formed  by  accumulation  of 
these  carriers  on  a  pixel  by  pixel  basis. 

Another  technique  involves  the  Charge  Injection  Devices 
<CID).  These  devices  are  surface  charge  devices  that 
collect  photon  generated  charges  and  store  them  in  MOS 
capacitors . 

Other  types  of  infrared  detectors  include  the  Lead 
Sulfide  (PhS)  Detectors,  Lead  Selenide  Detectors  (PbSe). 
Thermistor  Infrared  Detectors  and  Indium  Antimode  Charge 
Injection  Devices. 

One  tracking  technique  available  to  surveillance 
satellites  involves  the  use  of  the  doppler  effect.  The 
satellite  emits  a  signal  and  then  receives  the  signal '3 
return  after  it  reflects  off  a  given  target.  The  change  in 
frequency  as  a  result  of  the  relative  motion  between 
satellite  and  target  leads  to  a  calculation  of  the  targets 
velocity.  Ref  3  defines  the  return  frequency,  .  ,  detected 
by  the  satellite  as 

=  (1-u/c)/  l-(u/c)‘,  (2.29) 

where 

-  frequency  detected  if  both  were  at  rest, 
u  =  relative  separation  speeds, 
c  =  speed  of  light. 


Artificial  satellites  are  often  subject  to  a  variety  of 
disturbing  forces.  These  include  atmospheric  drag, 
variation  of  atmospheric  density,  solar  radiation  pressure, 
surface  charge  drag,  meteorite  impacts,  lunar  or  solar 
gravity  caused  perturbations  and  possible  encounters  with 
hostile  killer  satellites.  These  potential  disturbances  can 
act  seriously  to  impede  the  position  as  well  as  the  attitude 
stability  of  a  satellite.  Fortunately,  there  are  passive 
and  active  stabilization  systems  to  enable  recovery  from 


most  perturbations. 


III.  SATELLITE  COVERAGE  AND  COORDINATE  SYSTEMS 

This  chapter  includes  satellite  coverage,  spherical 
geometry  and  coordinate  systems.  A  coordinate  system  which 
enables  analytical  calculations  of  satellite  observer  and 
tracker  is  desired  for  at  least  a  single  case.  Several 
systems  are  explored  in  an  attempt  to  find  a  suitable  set  of 
reference  axes  to  establish  a  base  for  more  complicated 
analysis  by  computer.  A  simple  case  is  introduced  and  an 
appropriate  coordinate  system  is  presented. 

A.  GENERAL  SATELLITE  COVERAGE 

Basic  information  on  satellite  coverage  is  discussed 
here.  An  understanding  of  geometric  swath  width  and 
precession  is  required  for  a  detailed  study  of  satellite 
tracking . 

1 ■  Geometric  Swath  Width 

As  a  satellite  travels  around  the  earth  it  covers  a 
certain  amount  of  surface  area.  This  area  is  in  the  form  of 
a  spherical  cap  (figure  3.1).  As  the  satellite's  altitude 
above  the  earth  increases  the  area  of  the  spherical  cap 
increases.  The  area  of  the  spherical  cap  CRef.5)  can  be 
calculated  as 

A  =  2TTR1  (l-sin(90*-8)  )  <3.1> 
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where  R  is  the  radius  of  the  earth  and  9  is  the  angle  shown 
in  Figure  3.2. 

It  appears  from  this  that  a  high  altitude  orbit  is 
beneficial.  However,  it  must  be  considered  that  as 
satellite  altitude  increases,  satellite  sensor  accuracy 
decreases . 

If  a  target  is  at  point  'a'  <as  in  Figure  3.2)  it 
can  be  viewed  by  the  satellite  as  long  as  the  satellite  is 
between  points  1  and  2  on  its  orbit.  This  portion  of  the 
orbit  keeps  the  satellite  above  the  horizon  with  respect  to 
the  target.  The  amount  of  time  the  satellite  is  above  the 
horizon  is  calculated  CRef.  5]  to  be 

t  =  29 /w  (3.2) 

where  w  is  the  angular  velocity  of  the  satellite. 

As  the  satellite  progresses  in  its  orbit  the  spherical 
cap  of  coverage  moves  with  it.  As  it  moves,  the  cap  traces 
out  a  ribbon  around  the  earth.  The  width  of  the  ribbon 
(Figure  3.3)  is  the  geometric  swath  width  <GSW) . 

Figure  3.4  illustrates  the  GSW  as  the  very  dark  arc  on 
the  earth's  surface.  The  dotted  line  from  point  1  to 
point  2  is  tangent  to  the  earth  at  the  surface  location 
directly  beneath  the  satellite.  R  is  the  radius  of  the 
earth  and  h  is  the  satellite  altitude.  The  arc,  a(l>,  in 
Figure  3.5  is  one  half  the  dark  arc  in  Figure  3.4.  Simple 
geometry  proves  that  arc  a(l)  =  R9 .  Therefore  the 


GSW 


28R.  Figure  3.6  depicts  a  right  triangle 


From  the 


rules  of  trigonometry  cosine(Q)  =  R/(R+h).  Therefore, 

9  =  arccosine ( R/ C R+h ) ) 

It  follows  by  substitution  that 

GSW  =  2R < arccosine < R/ ( R *h )) )  (3.3) 

2 .  Ground  Track  and  Coverage 

As  a  satellite  travels  in  its  orbit  the- distance  it 
movgs  in  one  full  trip  around  the  earth  is 

D  =  2n ( R  +  h  >  (3.4> 

where  R  is  the  radius  of  the  earth  and  h  is  the  satellite's 
height  above  the  surface.  The  orbital  speed  is  calculated 
[Ref. 53  to  be 

-  4 . 2685 < R/ ( R *h > >  km /sec  (3.5) 

The  satellite  also  has  a  ground  speed  and  a  ground 
track.  The_  satellite's  position  can  be  projected  on  to  the 
earth's  surface.  This  point  on  the  earth's  surface  directly 
below  the  satellite  is  called  the  sub-satellite  point  (SSP). 
As  the  satellite  moves  its  SSP  forms  a  ground  track  which 
can  appear  as  simple  as  a  circle  about  the  earth  or  a  mere 
point  on  the  surface  or  they  can  appear  very  complicated. 

The  ground  speed  due  to  the  satellite's  motion  alone  is 

V  =  7 . 9053 (R/ (R+h) )km/sec  (3.6) 

However,  the  earth  does  rotate  at  about  0.267472 
km/sec  at  the  equator.  The  actual  ground  speed  depends  on 
the  inclination  of  the  orbit  plane  to  the  equatorial  plane. 
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For  all  cases  the  ground  speed  is  leaa  than  the  actual 
orbital  speed.  As  the  ground  track  is  formed  so  is  the 
ribbon  whose  width  is  the  GSW.  Because  of  the  earth 
rotation  beneath  the  satellite  the  swaths  may  overlap  at 
some  places  yet  never  cover  other  places.  If  the  earth's 
rotation  is  ignored  and  its  shape  is  assumed  spherical  a 
satellite  staying  in  one  plane  traces  the  same  circular 
swath  on  every  orbit. 

If  the  earth's  equatorial  bulge  is  taken  into 
account  the  plane  of  the  orbit  precesses  about  the  earth's 
axis  of  rotation.  This  precession  is  caused  by  the  non 
central  gravity  force  field  in  which  the  satellite  travels. 
It  is  measured  by  the  rate  of  drift  of  longitude  on  the 
surface  of  the  earth. 

For  an  elliptical  orbit  the  rate  of  precession  i3 

Si  =  9.95(R/a)  *  (cosine(i)/(l-ea)2  (3.7) 

where  SI  =  orbital  plane  precession  rate  in  degrees/day 
R  =  mean  radius  of  the  earth 
a  =  semimajor  axis 
e  =  eccentricity  of  ellipse 
i  =  angle  of  inclination  of  orbit 

For  the  special  case  of  a  circular  orbit  the 
precession  rate  is 

h  =  9 . 95 ( R/ ( R+h ) )  »•  cosine(i)  (3.8) 

If  the  effect  of  this  precession  is  ignored,  the 
motion  of  a  satellite  lies  in  a  plane  passing  through  the 


Kv? 


center  of  the  earth.  A  satellite  making  one  revolution  per 
day  at  0*  inclination  to  the  equator  is  stationary  over  a 
point  on  the  equator.  In  this  case  a  stationary  spherical 
cap  on  the  earth's  surface  is  the  only  area  covered. 
Commercial  communications  satellites  are  usually 
geostationary.  This  is  a  special  case  of  the  synchronous 
orbit  <24  hour  orbit) .  If  a  synchronous  orbit  is  inclined 
to  the  equator  the  ground  track  is  a  figure  eight. 

3 •  The  Spherical  Triangle 

At  this  time  all  the  simplifying  assumptions  are 
invoked  to  establish  a  base  for  further  analysis.  The 
satellite  is  traveling  in  a  circular  orbit  above  the 
equator.  Obviously,  transformations  to  other  orbits  are 
readily  available.  The  SSP  is  always  on  the  equator.  There 
is  a  target  traveling  on  a  steady  course  at  a  constant 
velocity  and  therefore  following  a  great  circle  heading. 

This  is  not  an  unreasonable  assumption.  It  is  normal  for  a 
ship  or  an  aircraft  to  keep  to  a  great  circle  heading  at  a 
most  efficient  cruising  speed  and  altitude.  The  target  is 
assumed  to  be  on  the  earth's  surface  or  very  near  to  the 
surface.  Figure  3.7  illustrates  the  paths  for  both  the 
surveilance  satellite  and  the  target.  The  target  track  is 
inclined  i*  to  the  satellite  track.  It  is  assumed  that  both 
can  be  considered  to  be  moving  on  the  surface  of  the  same 
sphere  with  transformations  readily  available  to  the  actual 


41 


:<■ 


satellite  position.  The  earth  is  assumed  to  be  a  perfect 
sphere  and  is  approximated  by  a  point  mass  at  the  center. 

The  two  paths  intersect  in  Figure  3.8  at  point  0. 

The  letter  S  represents  the  present  position  of  the 
satellite  and  the  letter  T  represents  the  present  position 
of  the  target.  These  three  points  <0,S  and  T>  on  the 
surface  of  the  sphere  define  a  spherical  triangle. 

Spherical  trigonometry  is  used  to  attempt  to  find  a  suitable 
coordinate  system  for  the  satellite  tracking  analysis. 

B.  COORDINATE  SYSTEMS 

Several  coordinate  systems  are  tested  at  this  time.  The 
set  of  reference  axes  is  a  relative  one.  It  is  assumed  that 
the  origin  of  the  coordinate  system  is  moving  with  the 
satellite  relative  to  the  earth.  Appendix  B  explains. 

1 •  Use  Variables  $  and  6 

Supposedly  any  two  independent  variables  from 
Figure  3.9  can  form  the  basis  for  a  state-space  model  for 
target  tracking.  The  variables  5  and  6  are  chosen  for  the 
first  attempt  to  form  a  suitable  coordinate  system.  The 
original  four  state  variables  x(l>  through  x(4>  are  defined 
below.  Derivatives  are  found  for  the  four  state  variables. 
In  an  attempt  to  establish  a  least-complicated  structure, 
however,  these  lead  to  other  state  variables  and  other 
derivatives:  i.e., 

x<l)  =  cosine(5> 


x  ( 2 ) 


sine<$) 


x<3)  9  cosine(6) 


x(4>  9  sine (6) 

Taking  derivatives  of  two  of  the  above  variables  yields 
x<l)  9  -x<2) 
xC3>  9  -x  C4) 

Using  the  trigonometric  law  of  sines  X(2)  can  be  defined  as 

x < 2 )  9  i  sine <t\) sine (t>  i  9  sine  (ip  sine  (6  ) 

[_  sine<6 )  J  sine(o) 

Since  in  this  case  a  and  sine(a)  are  constants  the  second 
identity  is  easier  to  use.  Sine(a)  is  replaced  by  the 
letter  A  and  sineC6>  is  replaced  by  x(4) .  This  yields 

X<2)  9  sine<T[>X(4>  /A 

The  same  procedure  is  used  for  X(4> .  The  law  of  sines 
provides  an  easy  alternate  form.  X(4)  can  be  defined  as 

x<4)  9  sine(T>A/sine(0> 

Table  3  summarizes  the  mathematics  of  expanding  the  state 
space.  It  is  soon  obvious  that  this  state  space  is  not 
readily  reducable  to  a  single  structure  by  increased 
dimenssion.  It  therefore  is  abandoned. 

2 .  Use  of  Variables  9  and  2 

Two  new  independent  variables  are  now  used  in  an 
attempt  to  generate  a  state-space  model.  As  before  the  sine 
and  cosine  functions  are  employed. 


X(l) 


sine<2) 


TABLE  3.  STATE  SPACE  FOR  VARIABLES  S  and  & 


XC1)  =  cosine (5)  X (2)  =  sine(S)  ! 

i 

X (3)  =  cosine ( & )  X(4)  =  s±ne(6)  j 

X<1>  =  -X(2)  and  XC4)  =  -X(3)  j 

X<2)  =  sin  sin <6  ) /sin <a>  *  sin <T^)  X (4)  /A 

X  (2)  =  cos<rpr\X(4> /A  ♦  sin  (TP  X  (  4  )  /  A 
Let  XC5)  =  cosO\>f\X(4>  and  X(6>  *  sin<q>X(4> 

(5)  =  -sin<T\)  CT\>  *X<4>  ♦  cos  (T\>T\X  <  4  >  ♦  cos<r\>  <T\>  (VpX(4> 
Let  X  (7)  *  sin <T\)  <T\>  1 X <4>  and  XC8)  =  cos(TpT\XC4> 
and  X  (9)  =  cos  (Tp  <fp  (T\>  X  <  4  > 

(7)  =  cos(f\)  CT\f  X<4>  ♦  sin(Tl)  (2T\)t\X(4)  ♦  sin (1\>  (TP  1  X < 4 )  I 
Let  X<10)  =  cos<T\)  <n/  X(4>  and  X(ll)  =  sin<rp  C2rpnX(4>  I 
and  X  <  12)  =  sinCfp  (r\)  *X(4>  1 


X(2) 


cosine(S) 
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X<1> 


sine ( 5  > 


X<2> 


cosine(l) 


XC3)  =  aine<8>  X(4)  =  cosineO) 

X<2)  =  -X(l)  X (4  >  =  -XC3) 

X(l)  =  sin  <T|)/  sin  (  8  )  X<3)  =  sin  («>  /  sin  < g)  \ 

(  sin(a)j  {^sin<T\)J 

XC1)  =  X < 3)  sin  (f\>  /  A  +  cos(TpT^X(3)  /  A 

Let  X  (5)  =  X  ( 3 >  s i n  (T\>  and  X<6>  =  cos  (T\)T\X  ( 3 ) 

X  <5)  *  X  ( 3)  sin (1\>  ♦  cos  <T\)1\X  ( 3) 

Let  X  <7)  =  XC3)sin(T\)  and  XC8)  =  cos  (T\>  QX  ( 3  > 

X ( 7 )  =  X <3) sin  (r\)  ♦  X C 3) cos Cf\)f\  =  X C 3 > cos <T\> t\;  <XC3>  = 

Let  X  <9 )  =  X  ( 3 > cos  (T\) r\ 

X  ( 9 )  =  X  <  3 )  s i n  <T\>  (T)_>  a  «■  X  ( 3)  cos  (T\> f\ 

Let  X(10)  =  X(3)sin(T\)  CO)  1  and  X(ll)  =  X(3)cos<n>r\ 

x<io)  =  x < 3 ) cos <r\>  a\>  *  X(3>sincri)  < 2i\) n_ 


ure  3.  10 


Projected  Track 
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Taking  time  derivatives  we  find 

p  »  . 

0  =  s  —  t  =  w<t>  -  y<s) 

To  obtain  a  state  space  interpretation  Y  and  X  and 
their  derivatives  are  defined.  Figure  3.13  shows  an  angle 
bisector  which  splits  9  and  ST'.  This  forms  two  right 
triangles.  Using  geometry  and  trigonometry  it  follows  that 

sine (9/2)  =  T'A/R 
AT'  =  ST'/2 
ST'  =  2Rsine<  9/2)  =  Y 
Y  =  2rcosine C 9/2) 9/2  =  Rcosine( 9/2) 9  = 

Rcosine (9-t) (w(t)-w(s) ) 

X  =  ST'  =  9R 
X  =  9R  ♦  R9 , 


but  R  =»  0; 


therefore. 


X  =  R9  =  R(w<t ) -w(s) ) 


Recall  R  is  a  constant  for  this  case  and  therefore  R 


is  zero. 


With  Y  and  X  defined  above  a  logical  state  space  is 
originated.  The  measurement  equation  is  also  important. 

In  order  to  define  a  measurement  equation 
Figure  3.14  is  constructed.  Note  that  points  0,S  and  T'  do 
not  form  a  straight  line. 

0T'  f  OS  +  ST' 

However,  in  terms  of  arc  lengths  the  following  applies; 


tiWTOWW;! 


o 


OT'  =  OS  +  ST' 

The  dotted  line  bisects  angle  a  and  therefore 

BT  =  TT'/2 

Using  geometry  and  trigonometry  it  can  be  shown  that 

TT'  =  20T ' sine (a/2 ) 


BT/OT  =  sine(a/2) 

From  Figure  3.10  it  follows  that 

OT'  =  2Rsine(  (p+t)/2) 

A  new  triangle  is  defined  in  Figure  3.15.  Using  simple 
geometry  it  is  obvious  that 

fi  =  w/2  -  a/2. 
and  by  simple  trigonometry 

DT  =  TT'sineCfi) 


DT '  =  TT'cosine(b) 


therefore 


tangent(€)  =  TD/SD 


SD  =  TD/ tangent (€ ) 


The  line  segment  ST'  is  solved  for 

ST'  =  SD  *  DT'  =  TD/ tangent <€ >  *  TT'cosine(fi) 

It  is  essential  to  express  ST'  in  terms  of  usable  variables 
ST'  =  TT '( s ine ( R ) /tangent (€ )  *  cosine(fl)) 
is  achieved  by  substitution  for  segment  length  SD . 

TT'  =  20T ' sine (a/2)  =  2 < 2Rsi ne ( ( p*t ) /2 ) ) s ine ( a /2 ) 
therefore 


ST'  =  4Raine  o+jJ  aine(a/2)  alne(fl)  ♦  coaine(B) 
v.  2  ,  tangent (€> 

Let  4Rsine( (a+T) /2)sine(a/2)  =  Q  for  convenience. 

Then  ST'  =  Q  l^sineCfl) /tangent  (€>  ♦  cosin Also, 

recall  that  fl  =  w/2  -  a/2 

ST'  =  Q  sine<Tr/2-q/2)  +  cosine<ir/2-a/2> 

\_"  tangent  ( € ) 

It  is  important  to  relate  X  and  Y. 

X  =  ST'  =  9R  and  Y  =  ST'  =  2Rsine<9/2> 

For  convenience  use  a  =  X/R 

Y  =  2Rsine(x/2R)  and  R  =  X/a  =  Y/2Rsine  (cr/2 ) 
therefore 

X  =  crY/2sine(cr/2> 

It  is  now  necessary  to  solve  for  €  which  indicates  bearing. 

Y  =  ST'  =  Qcosine < a/2 ) cotangent (€ )  Qsine(«/2) 

cotangent<€>  = _ Y _  _  _  -  sine(a/2) 

OcosineCa/2)  cosineva/2) 

Taking  the  arccotangent  of  each  side,  substituting  back  m 

for  Q  and  then  simplifying  leads  to 

€  =  arccotangent  sine^X_2R_) _  -  tangent  <  a/2) 

sine($/2)sine(a> 

This  is  the  measurement  equation.  The  final  result  is  a  two- 
variable  state  space  and  a  measurement  equation.  From  the 
measurement  equation  an  observer  can  be  modeled.  Table  5 


summarizes  the  results, 


State  Space 


arccot 


(3-14) 


IV.  OBSERVABILITY 


At  this  time  it  is  important  to  determine  observability 
of  the  target.  Assuming  the  satellite's  detectors  yield 
only  the  bearing  information  (as  developed  in  Chapter  3), 
it  is  necessary  to  obtain  the  target's  position  defined  with 
X  and  5.  Since  this  cannot  be  measured  directly  it  must  be 
estimated  in  some  way.  Analyzing  the  state  equations  from 
chapter  three  it  appears  that  if  ,^rand  are  given  the 
knowledge  of  X(0)  and  SCO)  (the  initialed  states)  are 
sufficient  to  determine  X(t>  and  S(t>  all  along  the 
trajectory. 

X(t)  *  X<0)  ♦  R  j  Ou^CiT)  -u^TO  ]d  c' 

'O  t 

S(t)  a  5(0)  +J  UJV  (  ir) d  L 

The  possibility  of  estimating  the  state  from  the  output 
measurement  under  the  assumption  that  the  system  dynamics  is 
completely  known  is  defined  as  observability.  Therefore 
the  observability  of  the  system  must  be  determined. 

All  changes  in  a  given  system  must  be  reflected  in  the 
output  for  the  system  to  be  observable.  Figure  4.1  shows  a 
simple  example.  The  states  X  and  X  are  estimated  from 
knowledge  of  the  output.  If  all  states  are  observable  for 
all  times  then  the  system  is  completely  observable.  This  is 
defined  in  detail  in  [Ref.  6) . 
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Figure  4.1  A  Simple  System 


The  test  for  observability  on  linear  systems  is  fairly 
straight  forward.  In  general  the  observability  of  linear 
systems  is  established  by  examining  the  linear 
independencies  of  the  columns  in  the  matrix  function 
C(t)5<t,tc),  where  C<t)  i3  the  measurement  matrix  and  $(t,t> 
is  the  state  transition  matrix.  For  the  linear  time 
invariant  system  this  leads  to  a  simple  rank  test  on  the 
appropriate  observability  matrix.  To  illustrate  this 
approach  we  may  consider  the  following  linear  system, 

X  =  AX 
Y  =  CX 

where  A  is  a  square  matrix,  X  is  an  n-dimensxonal  state 
vector  and  Y  is  an  m-dimensional  measurement  vector. 

Assuming  that  YCt)  is  differentiable  up  to  the  (n-l)th 
order,  we  have 

Y'=  CX  =  CAX 
Y" =  CAX  =  CA 1 X 


A  matrix  is  set  up  as  shown: 


where  Y  is  the  extended  measurement  vector  of  dimension  mn. 

The  matrix  Q  is  the  observability  matrix  and  has 
dimension  mnxn.  Clearly  if  rank  Q0  *  n,  X  can  be  uniquely 

A. 

solved  in  terms  of  the  measurement  Y.  It  should  be  noted 
that  differentiability  of  Y(t>  is  not  required  for 
observability  of  linear  systems  and  that  the  above  result 
can  be  derived  in  a  different  way.  The  presented  procedure, 
however  offers  the  possibility  of  extending  this  result  to 
the  observability  of  non-linear  systems;  namely  in  the 
equation,  Y  =  QttX,  any  subset  of  n  linearly  independent 


equations  can  be  considered  Y  =  QX,  0.  can  be  regarded  as  the 
Jacobian  J  of  Y  with  respect  to  X.  Hence  if  rank  J  =  n  the 
linear  time  invariant  system  is  observable. 

A.  OBSERVABILITY  OF  NON-LINEAR  SYSTEMS 

In  the  past  it  has  been  very  difficult  to  determine 
observability  of  non-linear  systems.  Recently,  however,  a 
number  of  new  methods  for  obsevability  testing  have  been 
developed . 

A  fairly  simple  test  that  determines  observability  in 
non-linear  systems  is  developed  in  CRef.  63.  This  new 
method  tests  observability  and  identifies  any  unobservable 
states  that  may  exist.  If  the  system  is  unobservable  than 
the  measurement  doesn't  provide  enough  information  for  state 


-  •  >  •  >  -  •  W  W.»V_*  ■'.A  u  *  *  *  '  -*  k  mVjk  **  ' 


estimation.  System  observability  and  state  estimation  are 
directly  related. 

A  non-linear  system  can  be  represented  as  follows: 

X(t>  =  f(X(t>,  u<t>,t>  (4.1) 

where  fO  is  an  n-function,  X  d  RA  ,  U  6.  if  . 

The  measurement  equation  is: 

Y  < t )  =  h(X(t> ,t>  (4.2) 

where  h()  is  an  m-function,  Y  £  R"\ 

Y(t)  must  be  differentiable  up  to  the  (n-1)  order  and 
U(t)  must  be  differentiable  up  to  the  (n-2)  order.  A  state 
X. (ta)  is  observable  at  t  if  knowledge  of  the  input  U(t) 
and  the  output  Y(t)  from  time  to  to  t  ,  enables  X  (te )  to  be 
determined . 

Equation  4.2  is  differentiated  to  provide  a  necessary 
system  of  non-linear  measurement  equations. 


Y  =  h(x) 

Y*  =  . 3  h  +  -)  h  J  x  =  hT  ♦  hAf  =  h,  (x.u) 
t  O  x  .)  t 

Y  3  J  h'  ♦  o  h,  Jt  +  -S  h ,  ~>u  = 

<3  t  OX-'X  J  u  ut 


hltr  *  h ,  *  f  ♦  h ,  -  hi(x,u,u') 

n ,  a 

Y  =  -  h,n-i>Ar  *  hIn.i,  .  .  . 


.  i*-ij  ,  (tv-2.) 

*  hc  »-ou  U- 


/  /)  -  2. 

) 


Define  an  mn  measurement  vector  Y  by 


and  an  an  function  H()  to  be 


H  <•)  = 


The  functional  relation  in  vector  form  is 

Y  =  H  <  x , v  >  (4.3) 

where  v<t>  is  a  function  of  ,  i  =  l,....(n-2>. 

With  respect  to  equation  4.3  the  question  of  observability 
can  be  developed  as  the  existence  of  an  inverse  of  this 
function.  Clearly,  this  is  related  to  the  Jacobians  of  H. 
However,  for  the  case  of  linear  systems  an  inverse,  if  it 
exists,  has  to  be  unique.  It  will  be  shown  late  that  here 
this  need  not  be  the  case. 

Two  conditions  must  be  met  for  this  non-linear  system  to 
be  observable  in  the  strict  sense.  These  two  conditions  are 
connectedness  and  univalence.  Connectivity  is  a  necessary 
condition  for  observability  in  the  strict  sense  and  a 
necessary  and  sufficient  condition  for  observability  in  the 
wide  sense.  Connectedness  is  satisfied  if  every  state  is 
connected  to  the  output  in  some  way.  Univalence  is 
satisfied  if  every  state  is  uniquely  determined  m  terms  of 
the  measured  output.  A  one-to-one  mapping  is  considered 
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univalent.  Both  conditions  together  are  necessary  and 
sufficient  to  establish  observability  in  the  strict  sense. 

To  explain  this  consider  Y(t>  as  expanded  in  a  Taylor 

series . 

Y  <  t )  =  y<tc>  +  y'(t;  )  <t-t  )  ♦  Y*y  <t,  )  (t-tc  >  1  ♦  .. 

+  <l/<n-l>  !  )y‘'ft'°  <tc  )  <t-tf  )  ^  *  r ( t ) 

Knowledge  of  the  measurement  trajectory  Y(t)  is 
equivalent  to  knowing  the  coefficients  and  remainder  in  this 
equation.  This  Taylor  series  expansion  is  considered  unique 
and  therefore  the  coefficients  are  also  unique.  The 
coefficients  are  the  elements  of  the  measurement  vector  Y. 
Therefore,  any  state  X,  <t>  is  observable  as  long  as  it  is 
connected  in  a  one-to-one  manner  to  an  element  of  Y.  If  the 
connection  is  not  one-to-one  the  system  is  observable  only 
in  the  wide  sense.  An  example  of  this  follows. 

X‘  =  Y 
X  = 

In  this  example  X  can  equal  either  the  positive  or 
negative  square  root  of  Y.  This  is  a  multiple  valued 
function.  Therefore  the  mapping  is  not  one-to-one.  Hence 
the  system  cannot  be  observable  is  the  strict  sense. 

As  mentioned  earlier  the  connectedness  condition  is 
satisfied  if  the  existence  of  an  inverse  of  the  function 
can  be  established.  According  to  [Ref.  6]  the  inverse 
function  is  considered  established  if  the  determinant  of  the 
nxn  Jacobian  of  H  does  not  equal  zero  for  all  x  and  v.  If 
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J  *  0  then  one  or  more  states  are  unobservable.  If  H  has 
more  than  one  right  inverse  then  uni valence  is  not  satisfied 
and  the  system  is  only  observable  in  the  wide  sense. 

However  this  will  suffice  for  this  case. 

B.  JACOBIANS  AND  NON-LINEAR  ANALYSIS 

Recall  from  chapter  three  that  two  possible  bearing-only 
measurements  have  been  found.  These  are  illustrated  in 
Tables  5  and  6.  Both  the  measurement  equations  (bearing- 
only)  are  highly  non-linear.  The  Jacobian  matrices  in 
Tables  7  and  S  are  defined  using  the  measurement  equations. 

These  Jacobians  are  necessary  for  observability 
analysis.  It  is  essential  to  ascertain  whether  or  not  their 
determinants  are  equal  to  zero.  The  Jacobians  in 
Tables  7  and  8  are  baaed  on  the  state  space  and  measurement 
information  available  in  Tables  5  and  6  respectively. 

The  calculations  for  rj  ,  <j  tj  /.j,  x ,  S and  O 
concerning  the  variable  X  are  in  Tables  9  through  11.  The 
result  of  the  calculations  for  V  ,  . j •')  /j ^  /~jJ  and 

<Ji) / .^  (f  concerning  the  variable  Y  proved  very  complicated  and 
are  therefore  not  used. 

It  is  obvious  that  the  determinants  of  these  Jacobians 
are  going  to  be  difficult  to  analyze.  Table  12  illustrates 
the  determinant  solution.  To  discover  whether  or  not  a  zero 
value  is  achieved  two  separate  methods  are  invoked.  First, 
a  simplification  is  utilized.  Recall  from  chapter  one  that 
a  satellite  is  restricted  regarding  the  amount  of  area  it 


can  cover.  Therefore  the  numerical  value  of  the  angles  $ 
and  t  are  restricted.  The  variable  X  is  composed  of  $  and  -r 
<X=R<8--t>>.  The  development  of  the  spherical  triangle  in 


TABLE  9 


JACOBIAN 


TABLE  11  JACOBIAN  -  X  -  PART  3 
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chapter  three  assumes  a  limit  to  angles  S,  r  and  a.  For  the 
purposes  of  this  work  the  range  of  possible  values  for  these 
angles  are  confined  as  follows; 

S:  0  (plus  epsilon)  radians  -  tt/4  radians 

t:  0  (plus  epsilon)  radians  -  tt/4  radians 

a:  0  (plus  epsilon)  radians  -  k/2  radians 

These  limitations  are  reasonable  for  thi3  particular 
case  and  they  enable  a  useful  simplification.  For  small 
angles  the  sine  of  that  angle  can  be  approximated  as  the 
angle  value  itself  (in  radians).  Table  13  verifies  the 
assumptions  made  here. 

TABLE  13.  SIMPLIFICATIONS 


w/8  =  S/2  at  S's  maximum 
w/8  =  .3926991  radians 
sin  ( ir/S)  =  .326834 


For  this  case  the  angle  a  is  held  constant.  Using  all 
the  above  assumptions  and  approximations  the  Jacobians  and 
their  determinants  are  greatly  reduced  in  complexity. 

Tables  14  and  15  summarizes  the  development  of  the 
simplified  Jacobian  and  determinant  for  the  measurement 
equation  containing  the  variable  X.  This  Jacobian  is  chosen 


because  it  is  the  least  complicated  of  the  two. 


TABLE  14 


SIMPLIFIED  VERSION 


The  second  method  involves  the  use  of  computer 
simulation.  The  determinant  (with  the  variable  X  as  before) 
is  calculated  using  no  assumptions.  The  angles  $  and  -r  are 
incremented  in  small  steps  from  0  radians  to  ir/4  radians. 

The  angle  a  (that  is  actually  constant)  has  been  given  five 
different  values  for  test  purposes.  The  value  of  the 
determinant  approaches  zero  but  never  actually  achieves 
zero . 

Note  from  table  15  that. 

•  Jj  *  _  _  a J  R*  _ _ _  (X$  -  §X)  <4.4> 

X1  -  4R5SX  +  4R 1 5 1 S 

If  X$  a  $X  the  determinant  of  the  Jacobian  is  zero  and 
the  system  is  not  observable.  From  Table  15  it  is  aDDarent 
that  the  system  is  not  observable  if  X  and  $  are  linearly 
dependent . 

The  various  calculations  and  computer  simulations 
indicate  that  the  system  initially  developed  in  chapter 
three  (Table  5)  is  observable  under  general  conditions. 
According  to  the  results  of  the  computer  simulation  the 
determinant  never  actually  reached  zero  for  the  particular 
test  values  (50  test  values)  chosen.  However  there  are 
certain  conditions  dependent  upon  the  initial  values  $  and 
t  when  the  system  is  not  observable.  Equation  4.1  defines 
these  circumstances.  Since  the  system  is  observable  for 
most  cases  it  is  possible  to  design  an  adaptive  observer. 
This  is  the  goal  of  Chapter  5. 


V. 


OBSERVERS 


At  this  point  in  the  satellite  tracking  project  it  is 
necessary  to  approach  the  state  estimation  problem.  In 
general,  a  process  can  be  characterized  by  the  state  and 
measurement  equations; 

X(t)  =  AX ( t )  (5.1) 

Z(t>  CX(t>  (5.2) 

where  X(t)  is  an  n-element  column  vector  representing  the 
states,  Z(t)  is  a  q-element  vector  of  measurements,  A  and  C 
are  nxn  and  qxn  matrices  respectively.  The  order  of  n  is 
greater  than  or  equal  to  the  order  of  q .  It  is  desired  to 
obtain  knowledge  of  the  states  of  X(t>.  However,  it  is  Z(t) 
that  is  measurable  and  therefore  known.  For  this  reason  it 
is  necessary  to  estimate  X(t).  The  device  for  estimating 
the  state  X(t>  is  the  observer. 

A.  LINEAR  LUENBERGER  OBSERVERS 

A  model  of  an  open  loop  observer  can  be  constructed  as 
follows ; 

Y ( t )  =  AY ( t )  (5.3) 

This  model  operates  as  does  the  process  or  plant.  The 
estimation  error  is 

e ( t )  =  Y ( t )  -  X(t)  (5.4) 
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This  is  differentiated  with  respect  to  time(t). 

e(t>  =  Y(t>  -  X(t) 

=  AY  <  t  >  -  AX ( t  > 

=  Ae(t)  (5.5) 

Therefore , 

e ( t )  =  e4'  e (0 )  =  0<t)e(O)  (5.6) 

If  the  eigenvalues  of  A  all  have  negative  real  parts 
then  e(t>  approaches  zero  as  t  approaches  infinity.  The 
error  signal  e(t)  decays  at  a  rate  determined  by  the 
location  of  the  eigenvalues  of  A.  The  rate  of  response  of 
the  process  is  also  determined  by  the  eigenvalues  of  A. 

To  make  the  decay  rate  independent  of  the  dynamic 
process,  often  a  state  model  is  'driven'  by  an  error  signal. 
In  this  case  the  observer  is  characterized  by; 

Y  ( t )  =  AY  ( t )  «■  G(Z(T)  -  CY  ( t  >  )  r5.7) 

and  the  error  equation  is 

e ( t )  =  Y ( t )  -  X(t) 

The  matrix  G  is  an  nxq  gain  matrix  which  can  be  selected 
to  determine  the  rate  at  which  y(t)  approaches  X(t>.  The 
observer  equation  can  be  written  as; 

Y  =  AY ( t )  ♦  G ( ( X ( t ) -CY ( t ) )  =  AY ( t )  ♦  GCe(t)  (5.8) 

The  driving  term  is  GCe(t)  and  its  purpose  is  to  drive 
the  estimate  Y  towards  X.  The  error  equation  is 
differentiated  with  respect  to  time. 

a(t)  =  Y ( t )  -  X(t) 


e(t> 


(A-GC) (Y(t)-X(t)) 


(5.9) 


e(t)  =  Fe  ( t  >  (5.10) 

where  F  =  A-GC 

e(t)  =  epT'  e ( 0 >  (5.11) 

As  expected,  the  results  depend  on  the  initial 
conditions.  The  gain  matrix  G  must  be  chosen  to  piace  tne 
eigenvalues  of  F  at  suitable  locations.  In  designina  an 
observer  it  is  usual  to  place  the  eigenvalues  of  F  where 
desired  and  then  determine  G. 

If  the  process  has  an  input  signal,  u,  then  the 
appropriate  equations  are; 

X(t)  =  AX(t)  *  BU(t)  <5.12) 

2 ( t )  =  CX(t) 

The  observer  for  this  system  is  characterized  by  the 
following  equation. 

Y ( t )  =  AY(t) +  BU(t ) *G(Z( t ) -CY ( t ) )  (5.13) 

The  error  response  is  the  same  as  in  the  above  and  is 
given  by  equation  5.11.  A  block  diagram  of  the  observer  is 
shown  in  Figure  5.1. 

When  designing  a  linear  or  non-linear  observer  it  is 
desirable  that  the  observer  error  become  small  rapidly  and 
that  the  observer  not  be  very  responsive  to  noise.  It  is 
very  difficult  to  meet  both  these  goals  since  they  seem  to 
conflict.  The  observer  must  also  be  supplied  with  a  set  of 
initial  conditions. 


If  some  of  the  states  of  a  given  system  can  be  measured 


there  is  no  reason  to  estimate  them.  An  observer  of  reduced 
dimension  can  be  designed  to  estimate  only  those  states  that 
cannot  be  directly  measured.  For  example,  if  there  are  n 
states  in  a  system  and  q  of  them  can  be  measured,  the 
observer  only  requires  <n-q>  states.  Figure  5.2  illustrates 
the  reduced  order  observer. 

B.  NON-LINEAR  OBSERVER 

For  the  design  of  the  non-linear  observer  the  techniaues 
developed  in  [Ref.  7]  can  be  applied. 

Declare  two  given  n-dimensional  vector  valued  functions 
of  time  to  be  Y(t)  and  8(t>.  Let  these  two  functions  be 
related  by  H<8)  as  Y<t>  =  H(8>.  The  vector  H(9>  is 
considered  known. 

Assume  that  each  ith  row  H(8>  is  differentiable  at  least 
once  on  all  parameters  8.  Therefore  an  nxn  Jacobian  matrix 
is  defined  as 


i 


I 


J  <  8 )  =  J  H(8 ) 

J  9 


(5.14; 


JL-  (  8  > 
J 


J  H  ^  (  9  ) 

"j' 


<  5 . 15  ) 


Y  ( t )  =  H  r<  8  ) 


*  J  H  <  8  >  J>  8  = 
6  8  J  c~ 


J  (  8  )  0  (  t  ) 


t5. lb  ) 
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Figure  5.3  ia  the  general  non-linear  observer  referred 
to  in  CRef .  71 .  The  transpose  of  the  Jacobian  is  represented 
as  JO).  A  string  of  parallel  integrators  are  denoted  as ( 
and  8^  (t)  represents  a  solution  of  YCt)  =  HC8)  for  any  given 
trajectory  Y  <t).  In  regards  to  this  figure  there  exists  a 

a 

positive  scalar  s  and  a  time  T>0  such  that  for  time  t>T 
<j|  8S  <t)  -  8,4  <  t  >  j  j  <  s)  if  the  following  three  conditions  are 
met . 

1.  K  is  positive  definite 

2.  The  magnitude  of  the  determinant  of  the  Jacobian  is 
bounded  both  from  above  and  away  from  zero  from 
below  for  the  particular  trajectory  Y,  (t) . 

3.  The  magnitude  of  Y^ (t)  ia  bounded  from  above. 

In  addition,  by  increasing  the  minimum  eigenvalue  of  K 
the  scalar  s  is  made  arbitrarily  small. 

An  estimation  error  signal  is  defined  as: 

e(t>  =  8  ;Ct>  -  8  ,Ct)  (5.17) 

Rearranging  provides; 

03<t>  =  e ( t )  +  8d(t) 

Referring  to  Figure  5.3  the  following  relation  is 
stated ; 

0,  =  -KJr(0s>  (HO;,)  -  H ( 0  ^  )  )  (5.18) 

For  notational  convenience  the  following  is  defined; 

Z(e,0d)  =  H(e*8  j  >  -  HOj)  =  HO,)  -  HO.)  (5.19) 

In  this  case  sigma  (Z)  does  not  refer  to  a  summation. 

It  is  simply  a  variable. 
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It  follows  from  the  above  calculations  that: 
e  =  03  -  %  -  -KJT(.9j)  Z(e,9d)  -  9j 


C5.20) 


The  tine-varying  Lyapunov  functions  of  the  error  signal 


are  defined  as: 


V(e, t>  =  *Z  <e,9j  JZte^g,*  > 

Since  the  determinant  of  the  Jacobian  is  bounded 


(5.21) 


ZCe^g^)  =  0  if  e  -  0  then  it  follows  that  V(e,t)  >  0  if 
e  =  0.  Further  calculations  are  summarized  in  Table  IS. 

The  minimum  eigenvalue  of  J<e  +  9.,)J  (e+0j)  is  denoted  by 
/)<,  which  is  always  greater  than  zero.  The  minimum 
eigenvalue  of  K  is  denoted  by  A  K  which  is  also  greater  than 


zero  so  both  eigenvalues  are  positive.  It  follows  that; 


V  <  -Af  /)K^Z<e,0d  >||  )Z(e,9d) 

Y,  *  7(9^)  9j 
9^  *  J  ( 9^  >''  Yj 


(5.24) 


(5.25) 


(5.26) 


Since  Y  and  J(0  )  are  both  bounded  J(0j)  and  8j  are 
also  bounded  for  all  t  >  0.  Therefore  \j  0j  i{  -  C(  and 
j  |  J  (  9j  )  |j  <  C  , . 

V  <  -  Af,\|Z(«,0i  >|j*  +  C,C,||Z(e,ej  )|j 
V  <  -(wm7  |Z(e,9j  );|  -C  /„  2,-\v.W~)  *  +C,  C,/ 2  ,U  l  vr  (5.27) 
C.  DEVELOPMENT  OF  SATELLITE  OBSERVER 

For  the  design  of  the  satellite  observer  (the  system 
initially  presented  in  Table  5)  the  state  6  as  described 


above  now  represents  both  X  and  5  as  equation  5.28  shows: 


9  =  X, 


(5.28) 
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TABLE  16 


LYAPUNOV  FUNCTION  OF  THE  ERROR  SIGNAL 


Mil'll  & 

T~ 


06  - 


1  (e.  -<?J )  J  (e-Cj  )\sMe^)2  le,  &J 


V  =  -Zr(e,Pd\  )!>/,,  ]-<sJ  Jhj)Z(e,Sn) 


The  results  developed  earlier  in  this  chapter 
(part  B)  are  now  applied  to  the  particular  problem  of  the 
satellite  tracker.  The  Jacobian  of  Table  7  is  redefined  in 
equation  5.31. 

It  is  intuitively  obvious  that; 

Y  =  H  ( X ,  5  ,  U  >  =  H  <  8  ) 

The  further  development  of  both  H(8>  and  8  is  presented 
in  table  17.  The  information  in  Table  17  leads  to  the 
design  of  an  observer  for  the  satellite.  This  design  is 
depicted  in  Figure  5.3.  Figure  5.4  shows  the  model  used  to 
check  observer  results.  Both  the  observer  and  model  are 
simulated  using  the  FORTRAN  programming  language.  The 
results  are  discussed  in  the  next  chapter. 


TABLE  17  CALCULATIONS  FOR  OBSERVER  DESIGN 
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VI.  SIMULATION  RESULTS 

The  basic  observer  designed  in  Chapter  5  is  simulated  to 
prove  authenticity.  Recall  that  the  original  state  and 
measurement  equations  were  developed  in  Chapter  3.  Due  to 
the  satellite's  limited  GSW,  several  parameters  (including 
angles  $  and  t )  are  confined  in  their  maximum  attainable 
values.  Therefore,  these  parameters  were  kept  within 
specific  ranges.  The  angle  alpha  (a)  is  taken  to  be  a 
constant  30*.  This  value  is  selected  for  convenience. 

The  angular  velocity  of  the  target  ('....)  is  chosen  to  be 
0.0005  rad/sec.  This  sets  the  target  speed  at  about 
2  mi/sec.  The  angular  velocity  of  the  satellite  is 

chosen  to  be  0.00045  rad. sec  which  yields  a  speed  of  about 
1.8  mi/sec.  These  are  arbitrary  but  not  unrealistic  values. 
A.  THE  BASIC  CONTINUOUS  OBSERVER 

Figure  6.1  depicts  the  basic  simulation  block  diagram. 
The  initial  conditions  5  and  X  are  0.2  radians  and  1000km 
(621.5  mi)  respectively.  These  values  are  arbitrary  but 
within  the  acceptable  range.  Recall  that  X  =  R(3S-t). 

Figure  6.2  is  a  block  diagram  depicting  the  basic 
observer  design  that  is  simulated.  Note  that  Figure  6.2  is 
virtually  identical  to  Figure  5.5. 


►  y 


The  equations  representing  the  various  elements  m 
Figure  6.2  were  initially  developed  in  Chapters  3,4  and  5. 
Table  IS  summarizes  these  equations. 

The  gain  matrix  (as  developed  in  Table  17)  is  chosen  to 
be  a  diagonal  matrix  of  the  form: 

K  =  k,  0 

!  I 

j  0  kij 

The  elements  k  and  k  both  must  be  greater  than  zero. 

A  diagonal  matrix  was  chosen  because  it  is  fairly  easy  to 
work  with  and  it  allows  the  error  changes  to  be  monitored. 

There  are  two  ways  to  describe  the  observer  to  be 
simulated.  It  can  be  described  in  the  continuous  domain  or 
in  the  discrete  domain. 

The  observer  was  initially  represented  as  a  continuous 
model.  A  forth-order  Runge-Kutta  method  is  used  for 
integration.  The  continuous  model  was  used  in  order  to 
eliminate  any  possible  errors  that  could  be  introduced  as  a 
result  of  discretization.  However,  simulations  for  both 
versions  produce  no  real  differences. 

B.  THE  DISCRETE  OBSERVER 

With  the  advent  of  micro-computers  a  discrete  observer 
acting  as  a  one-step  predictor  is  a  more  realistic  choice. 
In  this  model  the  foreward  difference  approximation  is 
employed . 

The  basic  definition  is; 

XCt  )  =  X<t.  )  *  1 1  X<t.  ) 


85 
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Full  Non-Linear  Observer 


First  the  original,  highly  non-linear  observer 
(figure  5.5)  was  simulated  using  a  FORTRAN  program  on  the 
IBM  370  system.  All  attempts  to  find  values  of  gains  and 
integration  step  for  which  observed  states  converge  were 
unsuccessful.  Typically  the  observer  enters  steady  state 
with  very  small  errors  e  and  e^;  however,  the  estimated 
states  were  far  from  actual  values.  It  is  possible  that 
this  behavior  is  the  consequence  of  the  sinusoidal  periodic 
functions  (sin(x/2R)  and  sin($/2>)  in  7  and  .  Recall  from 
Chapter  4  that  this  observer  is  not  shown  to  be  observable 
in  the  strict  sense.  Univalence  may  not  exist  due  to  the 
reiterative  nature  of  the  sine  function.  In  order  to  avoid 
periodicity  an  approximation  of  a  non-linear  observer  is 
considered . 

The  simplified  observer  has  the  same  form  as  the 
full  observer  (see  Table  18),  but  the  measurement  equations 
and  the  corresponding  Jacobian  are  approximated  as  developed 
in  Chapter  4  (and  repeated  in  Table  19) . 

The  simplified  observer  simulated  was  being  very 
insensitive  to  changes  in  X  and  $.  The  observer  stabilized 

A  A 

with  a  small  error  (e ,  and  e j>  while  X  and  $  remained 
significantly  different  from  X  and  5  respectively.  The 
exact  reason  for  this  behavior  is  unknown.  However  it  is 
theorized  that  the  relative  smoothness  of  the  arccotangent 
curve  (arccot  ;i)  is  responsible  for  this  behavior. 
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TABLE  19 


SIMPLIFIED  OBSERVER 


Measurement  equations: 

set  s  =  sin1 <a/2>  and  a  =  sin(a) 
h  (X,$)  =  arccot  x  -  2R$S  ' 


V  R$a 


h  ;  <X,5)  =  -aR, 


X$  -  IX 


l  X'  -  4R5SX  *  4R1$1S 


Jacobian : 

J,>  =-v/.X 


-aR 


X'  -  4R5SX  ♦  4R 2  $ 1 S 


=  •  /  5  = 


aR 


/.'X  = 


X1  -  4R3>SX  +  4R1$1S 
aR 


X'  -  4R5SX  +  4R1$1S 
2X$ ( X-2R5S  >  ♦  S(4R'$,S-X1) 


~  /  A.,S 


X1  - "4RSSX  ♦  4R'?'S 
aR 


X‘  -  4R1SX  +  4R 1  $ 1  S 
X(4R1$‘S-X1)  -  4ESX$(X-2R$) 


X1  -  4R5SX  ♦  4R1$15 
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In  order  to  circumvent  this  problem  the  measurement 
signal  is  redefined. 

2 .  Redefined  Non-Linear  Model 

A  new  measurement  variable  is  defined  as: 

<t )  =  x(t)/i(t)  =  Rsin  <  o )  (cot  ( ■'>  ( t )  )  +  tan  (a/2)) 

If  the  value  of  a  is  known  and  y<t>  is  measured: 
then  u(t)  can  be  calculated.  Taking  into  account  the 
assumed  range  for  §(t)  and  X(t>  and  the  values  for 
parameters  R  and  a,  the  expression  X/ ( RSsin (a) ) -tan ( a/2 ) 
never  achieves  a  very  large  value.  Therefore  (t)  is 
bounded  away  from  zero.  This  implies  that  cotanC  (t))  has  a 
bounded  value. 

The  derivative  of  ,->-Ct)  is: 

£-(t)  =  Rsin(a)(-"(t>/sin'(;(t))) 

where  r"<t)  and  '  (t)  are  measured  values  and  (t)  is 

‘  / 

evaluated . 

The  redefined  non-linear  model  and  observer  are 
presented  in  Figure  6.3.  It  is  noted  that  this  observer  has 
the  same  form  as  the  observer  given  on  the  block  diagram  m 
Figure  6.2.  Observer  equations  are  shown  in  Table  20. 

All  attempts  to  find  the  power  K, ,  K  .  and  ..t  (the 
time  increment)  that  would  force  the  states  to  converge 
failed.  A  satisfactory  theoretical  explanation  for  this 
behavior  was  not  established.  However,  one  possible 
explanation  is  that  e  and  e  are  extremely  sensitive  to 


TABLE  20 


MEASUREMENT  EQUATIONS  -  REDEFINED  NON-LINEAR  MODEL 


Measurement  equations'. 

h  i  <  X  ,  S )  =/UX,i)  =  X 

S 

hv(X,35)  =  ^--(X.S)  =  XS  -  SX  =  1  ,X  -  SX 

s  i.  s 

Jacobian: 

J ,  i  =  _^_h_  =  ;  J,.  =  h  •  =  -X_.  =  -1  X 

X  "s’"  S  S'  Si 

J:,  =  -3_h =  -1'  S 

X  S' 

J  =  .  h  =X( S  > '  -  2$ ( XS-SX )  =  -i  :  XS  *  2SX 

S  S  ‘  ‘  S’ 


=  1 _ ,  2SX 

S'  s 


X  | 


Observabi 1 i ty ; 


jjj  =  J,,  J,,  -  Jlt  J.,  =  X-  2SX  -  X  ;  -  1_.  SX  =  1  .  SX 

$  *  ^  f  ^  ^  ^  ^ 


J!  =  O  if  IX  -  X  =  O 


X  =  $ 


X  /  $ 


i  Note  that  these  conditions  were  defined  previously. 


Analysis  of  the  model  with  redefined  measurements 
reveals  that  the  form  of  the  equations  allow  redefinition  of 
state  variables.  This  may  reduce  the  sensitivity  of  the 
model . 

A  new  state  vector  CU,V]'r  is  defined  as: 

U  =  X/$  and  V  =  1/S 

Table  21  illustrates  further  development. 

The  block  diagram  in  figure  6.4  shows  the  discrete 
observer .  A  simulation  of  the  discrete  observer  was  run 
with  =  0.0005  and  =0 . 00045 .  Several  different  gains 
were  used  to  test  convergence.  Convergence  is  obtained  for 
the  following  ranges: 

0  <  K,  At  <  2  and  0  <  K  •  At/ j  J  _  .  j  <  2 

A  /*» 

Arbitrary  initial  values  for  U  and  V  were  employed. 
It  seemed  natural  to  assume  Un  =  0 ,  X  =0  and  V  =  1  (5  =1 
is  within  the  allowable  value  for  $). 

To  examine  the  observer's  ability  to  track  varying 
angular  speeds  it  is  assumed  that  ii'rand  ui  are  sine  waves. 

=  A  sin(Dt) 

=  B  sin(Dt) 

Different  amplitudes  and  periods  are  tested  and 
successful  tracking  is  achieved.  Appendix  C  contains  the 
resulting  computer  graphs. 


TABLE  21 


OBSERVABILITY  CONDITIONS  AND  DISCRETIZED  OBSERVER 


Measurement  equations: 


/  U  U  ,  V  >  =  h  ,  <  U  ,  V  >  =  U 
?_<-(U,V)  =  h  ,_  ( U  ,  V )  =  V(X  - 

Jacobian : 


J , ,  =  -Jh,  =  1 
o  U 

:  J,;  *  lJl  = 

«■  V 

0 

j_ ,  s  '  h:  =  -SV  ; 

■  u 

J,,  =  vh.  =  X 

V<i 

- 

$u 

Observability  conditions 

are  the  same: 

• J !  =  X  -3U  P  0 

|  j 

X  -  $X  *  0  X 

r 

$ 

/’ 

/• 

Continuous  observer: 

U  =  K,  J  , 

e  i  *  KJ  e  . 

II 

'> 

e  ,  ♦  K  J;.  e  . 

The  discretized  observer 

becomes : 

e ,  K  =  i  (Ut  ,Vr  >  - 

h,  (U-  , V,  ) 

u , 

e  ?  v  =  /  (U  .  Vr  )  -  h  ,  (U 

.  .V, )  =  •  .  -  vf (X 

- 

$  U,  ) 

U ,  ,  ,  =  U,  ♦  K,  ite„ 

♦  K  .  t  •  J  •_  i  (  U 

.V , 

,  >  e  . 

Vm,  =  V.,  * 

Ok  ,Vf  )e, 

As  expected  the  amplitude  of  the  sine  wave 
representing  state  <A/D;B/D>  which  can  be  tracked  is 
inversely  related  to  the  period  of  oscillation. 

4.  Model  Decouolinq 

Analysis  of  the  error  equations  provide; 

e.  =  ;*<U,V>  -  h,<U,V>  =  U-U 

A.  «  ,  \  i 

e,«  ,WU,V)  -  h  j.  <  U  ,  V  )  =  V(X-$U) -V(X-$U) 

=  X<V-V) -$(UV-UV) 

The  block  diagram  in  Figure  6.4  indicates  the 
possibility  of  decoupling  the  observer  into  two  first  order 
models  by  letting  J^,  equal  zero  (J,>  is  already  zero)  . 

Figure  6.5  shows  a  block  diagram  of  the  discrete 
subsystem.  The  discrete  system  of  Figure  6.6  is  certainly 
stable  if  it's  eigenvalue  lies  within  the  unit  circle. 

i  1  -  a  i  <  1  ^  O  <  a  <  2  0  <  K  it  <  2 

I  ! 

Under  this  condition  the  system  has  a  steady-state 
error  of  zero  for  a  constant  input  u. 

Hence  the  system  acts  as  positioning  servo-mechan  is 
and  is  capable  of  tracking  the  relatively  slow  varying 
inputs  with  very  small  errors. 

5.  Subsystem _ Two. 

Assuming  that  subsystem  1  is  tracking  u  <  u<t>  = 
u<t))  the  error  e  c  can  be  approximated  as: 

,  .v 

e  *  *  (X  -  $U>  (V  -  V) 

The  observer  is  reduced  to  the  gradient  algorithm 
where  the  gradient  is: 


Figure  6.7  illustrates  the  observer.  This  result  is 
also  in  accordance  with  the  method  proposed  in  CRef .  6] . 

*  *  'V. 

Defining  <X-SU >V  as  a  new  measurement  the 

»  .  A 

corresponding  function  h  is  h=<X-$U)V.  The  Jacobian  is 
J  =  X-5U. 

In  order  to  simplify  the  stability  analysis  the 
Jacobian.  J,  is  replaced  with  the  signum  function, 
signum(J).  This  is  sufficient  to  guarantee  the  negative 
feedback  of  subsystem  2.  Subsystem  2  is  analogous  to 
subsystem  1  and  will  be  stable  if  0  <  K .  it  <  2. 
Consequently  this  system  also  performs  as  positioning  servo¬ 
mechanism  and  will  track  relatively  slow  varying  inputs  V(t> 
with  very  small  error. 

The  presented  analysis  leads  to  a  definition  of  the 
decoupled  observer.  This  observer  is  depicted  in  Table  22. 
Appendix  D  contains  the  computer  simulation  results. 

Figure  6.8  is  essentially  the  same  as  Figure  6.4 
with  J^<  =  0  and  JLi  =  signum(J2i  ).  Therefore  it  is 

possible  to  consider  this  algorithm  a  simplification  of  the 
algorithm  developed  in  CRef.  71.  However,  it  should  be 
pointed  out  that  the  decoupled  observer  can  be  derived 
directly  through  analysis  of  the  two  linear  subsystems 
presented  above. 


96 


In  the  decoupled  model  it  is  easier  to  control  gains 


since  the  JLi  element  does  not  affect  stability.  It  is 
noted  that  slightly  more  oscillatory  behavior  is  present 
at  the  beginning  of  the  simulation.  This  is  not  unusual  for 
gradient  type  algorithms. 

Satisfactory  results  for  the  product  Kit  =  0.1  are 
expected  if  the  observer  is  analyzed  from  the  sampling 
theorem  point  of  view.  This  presupposes  that  both 
continuous  subsystems  posses  only  one  time  constant  =  1/K. 
The  sampling  theorem  requires  that  t  <  t~/2  and  4tK  <  0.5. 

Experience  has  shown  that  the  best  results  for  the 
one-step  integration  are  obtained  if  the  integration  step 
(sampling  time)  is  ten  times  less  that  the  smaller  time 
constant  in  the  system.  This  means  4.t  *  0.1  (1/K)  and 
*tK  *  0.1. 

The  sampling  theorem  is  violated  if  AtK=l.  As  a 
consequence  the  observer  requires  much  more  time  to  extract 
sufficient  information  from  the  measurement  and  to  start 
state  tracking.  The  simulation  results  demonstrate  that 
when  the  product  atK  equals  unity  convergence  starts  after 
fifty  seconds.  This  is  approximately  ten  times  slower  from 
previous  experiments  where  the  sampling  theorem  is 


V  j 
<:> 
/„  J 
v/; 
>  \ 


satisfied . 

In  the  case  of  itK  =  0,01  the  system  pole  is  very 
close  to  the  unit  circle.  Therefore,  convergence  is 


As  demonstrated  by  the  various  computer  generated 
graphs,  the  simulation  results  for  this  system  are 
satisfactory.  Therefore  the  selected  coordinate  system  is 


proven  acceptable 


VI.  SIMULATION  RESULTS 


The  basic  observer  designed  in  Chapter  5  is  simulated  to 
prove  authenticity.  Recall  that  the  original  state  and 
measurement  equations  were  developed  in  Chapter  3.  Due  to 
the  satellite's  limited  GSW,  several  parameters  (including 
angles  S  and  -r>  are  confined  in  their  maximum  attainable 
values.  Therefore,  these  parameters  were  kept  within 
specific  ranges.  The  angle  alpha  (a>  is  taken  to  be  a 
constant  30*.  This  value  is  selected  for  convenience. 

The  angular  velocity  of  the  target  (a-*.)  is  chosen  to  be 
0 . 0005  rad/sec.  This  sets  the  target  speed  at  about 
2  mi/sec.  The  angular  velocity  (u~>  of  the  satellite  is 
chosen  to  be  0.00045  rad. sec  which  yields  a  speed  of  about 
1.8  mi/sec.  These  are  arbitrary  but  not  unrealistic  values. 
A.  THE  BASIC  CONTINUOUS  OBSERVER 

Figure  6.1  depicts  the  basic  simulation  block  diagram. 
The  initial  conditions  S  and  X  are  0.2  radians  and  1000km 
(621.5  mi)  respectively.  These  values  are  arbitrary  but 
within  the  acceptable  range.  Recall  that  X  =  Rd-T). 

Figure  6.2  is  a  block  diagram  depicting  the  basic 
observer  design  that  is  simulated.  Note  that  Figure  6.2  is 
virtually  identical  to  Figure  5.5. 


F i gure  6.  1 


Basic  Simulation  Block  Diagram 


The  equation*  representing  the  various  eleaents  in 
Figure  6.2  were  initially  developed  in  Chapters  3.4  and  5. 
Table  18  suaaarizes  these  equations. 

The  gain  Matrix  (as  developed  in  Table  17)  is  chosen  to 
be  a  diagonal  matrix  of  the  form: 

k  =  r  k,  oi 


The  elements  k  and  k  both  must  be  greater  than  zero. 

A  diagonal  matrix  was  chosen  because  it  is  fairly  easy  to 
work  with  and  it  allows  the  error  changes  to  be  monitored. 

There  are  two  ways  to  describe  the  observer  to  be 
simulated.  It  can  be  described  in  the  continuous  domain  or 
in  the  discrete  domain. 

The  observer  was  initially  represented  as  a  continuous 
model.  A  forth-order  Runge-Kutta  method  is  used  for 
integration.  The  continuous  model  was  used  in  order  to 
eliminate  any  possible  errors  that  could  be  introduced  as  a 
result  of  discretization.  However,  simulations  for  both 
versions  produce  no  real  differences. 

B.  THE  DISCRETE  OBSERVER 

With  the  advent  of  micro-computers  a  discrete  observer 
acting  as  a  one-step  predictor  is  a  more  realistic  choice. 
In  this  model  the  foreward  difference  approximation  is 
employed . 

The  basic  definition  is; 


TABLE  18  OBSERVER  EQUATIONS 
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Full  Non-Linear  Observer 

First  the  original,  highly  non-linear  observer 
(figure  5.5)  was  simulated  using  a  FORTRAN  program  on  the 
IBM  370  system.  All  attempts  to  find  values  of  gains  and 
integration  step  for  which  observed  states  converge  were 
unsuccessful.  Typically  the  observer  enters  steady  state 
with  very  small  errors  e,  and  ee;  however,  the  estimated 
states  were  far  from  actual  values.  It  is  possible  that 
this  behavior  is  the  consequence  of  the  sinusoidal  periodic 
functions  <sin(x/2R>  and  sin($/2))  in  7?  and  '1  .  Recall  from 
Chapter  4  that  this  observer  is  not  shown  to  be  observable 
in  the  strict  sense.  Uni valence  may  not  exist  due  to  the 
reiterative  nature  of  the  sine  function.  In  order  to  avoid 
periodicity  an  approximation  of  a  non-linear  observer  is 
considered . 

The  simplified  observer  has  the  same  form  as  the 
full  observer  (see  Table  18),  but  the  measurement  equations 
and  the  corresponding  Jacobian  are  approximated  as  developed 
in  Chapter  4  (and  repeated  in  Table  19) . 

The  simplified  observer  simulated  was  being  very 
insensitive  to  changes  in  X  and  $.  The  observer  stabilized 

A  A. 

with  a  small  error  (e,  and  e^)  while  X  and  $  remained 
significantly  different  from  X  and  $  respectively.  The 
exact  reason  for  this  behavior  is  unknown.  However  it  is 
theorized  that  the  relative  smoothness  of  the  arccotangent 
curve  (arccot  ?7>  is  responsible  for  this  behavior. 


TABLE  19  SIMPLIFIED  OBSERVER 


Measurement,  equations: 

set  s  =  ain1 (a/2)  and  a  =  sin(a) 

h  <X ,  5)  =  arccot  ( x  -  2R5S  ^ 

l  RSa  J 

h;(X,5)  =  -oRL _ X5  -  fX _ 

i  X*  -  4R5SX  +  4R151S/ 

Jacobian : 


Jm 

=  Jy/Jx 

= 

-aR 

5 

( 

X1  - 

4R5SX  - 

4R*5* 

S 

=  «J  y  /J  s 

3 

aR 

X 

X* 

-  4R5SX  ♦ 

4Ra  5 

aS 

J  2.  \ 

*  'J'y/J  x 

— 

aR 

* 

X* 

-  4R5SX 

♦  4R* 

5aS 

2X5 <X- 

2R5S)  +  5 

i(4R*5 

*s-xa  > 

X1  - 

4R5SX  * 

4Ra5* 

S 

=  -JV/S* 

= 

aR 

* 

X1  -  4R5SX  +  4R'5aS 


X(4R‘$iS-X>)  *  4ESX5 ( X - 2R5 ) 
X'  -  4R5SX  >■  4R'5a5 


In  order  to  circumvent  this  problem  the  measurement 


signal  is  redefined. 

2.  Redefined  Non-Linear  Model 

A  new  measurement  variable  is  defined  as: 

/-'-(t)  =  x(t)/2(t)  =  Rsin  (a)  (cot  (  ^<  t )  )  +  tan(a/2)> 

If  the  value  of  a  is  known  and  ^/(t)  is  measured: 
then  u(t)  can  be  calculated.  Taking  into  account  the 
assumed  range  for  £(t)  and  X(t>  and  the  values  for 
parameters  R  and  a,  the  expression  X/ (RSsin (a ) > - tan (a/2 ) 
never  achieves  a  very  large  value.  Therefore  ^ (t)  is 
bounded  away  from  zero.  This  implies  that  cotan(^(t) )  has  a 
bounded  value. 

The  derivative  is: 

.O(t)  *  Rsin(a)  (-^(t>/sinJ  (?(t>  )  ) 
where  ^(t)  and  >?( t)  are  measured  values  and  ,u(t)  is 
evaluated . 

The  redefined  non-linear  model  and  observer  are 
presented  in  Figure  6.3.  It  is  noted  that  this  observer  has 
the  same  form  as  the  observer  given  on  the  block  diagram  in 
Figure  6.2.  Observer  equations  are  shown  in  Table  20. 

All  attempts  to  find  the  power  K(  ,  K  .  and  at  (the 
time  increment)  that  would  force  the  states  to  converge 
failed.  A  satisfactory  theoretical  explanation  for  this 
behavior  was  not  established.  However,  one  possible 
explanation  i3  that  e  and  e  are  extremely  sensitive  to 


TABLE  20 


MEASUREMENT  EQUATIONS  -  REDEFINED  NON-LINEAR  MODEL 


Measurement  equations: 


h,  < X ,  a)  =  /^<X,a>  =  J L_ 


h2cx,5)  =  M<x,a>  =  xa  -  ax  =  1  [x  -  axj 

a  *  at  a ! 

Jacobian: 


Jh  =  -J  h  .  =  __1__  ;  Jlt  =  J  h,  =  _rX_  =  -JL_  X_ 

^jx  a  a  aj  aa 


Jn  =  <-J  h  t 

0  X 


•it  gj 

a1 


j4l  =  -j  h.  =  x<a>j  -  za<xa-ax) 


•  zl-C 


xa  +  2ax 


=  i__/2ax_-  x 

(  a 


Observability ; 


j  J  j  =  J, i 


j{i_  ji(  =  i  2 ax  -  x 

a1 1  a 


i_  ax  wax. 

a3  a  a‘i  a 


i  j  j  =  o  if  ax.  -  x  =  o  x  =  a, 
a  x 


Note  that  these  conditions  were  defined  previously, 


change  in  states  which  results  in  an  extremely  slow 


convergence  with  a  highly  oscillatory  transient  response. 


Analysis  of  the  model  with  redefined  measurements 
reveals  that  the  fora  of  the  equations  allow  redefinition  of 
state  variables.  This  may  reduce  the  sensitivity  of  the 


model . 


A  new  state  vector  CU,V3  is  defined  as: 


U  =  X/5  and  V  =  1/5 

Table  21  illustrates  further  development. 

The  block  diagram  in  figure  6.4  shows  the  discrete 
observer.  A  simulation  of  the  discrete  observer  was  run 
withes  0.0005  and  *0.00045.  Several  different  gains 
were  used  to  test  convergence.  Convergence  is  obtained  for 
the  following  ranges: 

0  <  K,  4t  <  2  and  0  <  KtAt/JJ^j  <  2 

A  A 

Arbitrary  initial  values  for  U  and  V  were  employed. 

A  A 

It  seemed  natural  to  assume  U«  =  0r  X  =0  and  V  =  1  (5  =1 
is  within  the  allowable  value  for  5) . 

To  examine  the  observer's  ability  to  track  varying 
angular  speeds  it  is  assumed  that  U^.and  w*  are  sine  waves. 

=  A  sinCDt) 

C>Js  -  B  sin(Dt) 

Different  amplitudes  and  periods  are  tested  and 
successful  tracking  is  achieved.  Appendix  C  contains  the 
resulting  computer  graphs. 


•  *->  v  -j.  ■  ■  7’.-  v  v  -\ 


Redefined  Model 


TABLE  21 

OBSERVABILITY  CONDITIONS  AND  DISCRETIZED  OBSERVER 

Measurement  equations: 

/  U  U ,  V  >  =  h ,  <  U , V)  =  U 
Xi<U,V>  =  (U, V)  =  V<X  -  1U) 

Jacobian : 

I 

J»,  *  *  1  :  J,2  -£h.  =  o  | 

JU  OV 

I 

J».  =  i  ht  =  -SV  ;  JiZ  =  rJJVj  =  X  -  5U  i 

J  U  JV 

I 

I 

Observability  conditions  are  the  same:  j 

| 

I J I  =  X  -SU  AO  X-SX/O  <5 

^  *  <A'  I 


Continuous  observer: 

U  =  K,  J)(  e,  -  K.J,,  e, 

V  =  KjJu  e,  K  iJu  &i 
The  discretized  observer  becomes; 

mlK  ~  ;U<U«  ,Vk>  -  h,  (Up;  ,  V^  >  =/-*.-  Uc 

e^m-\)  -  ht<U„,v\>  =  Vr<X  -  S 

U„fl  *  UK  ♦  K |  <i telfc  ♦  K  ,  at  •  JiltUMV,)! 
VK,,  =  V„  ♦  K.itJ^  <U*  ,VK)eiCt 


c> 


As  expected  the  amplitude  of  the  sine  wave 
representing  state  <A/D;B/D>  which  can  be  tracked  is 
inversely  related  to  the  period  of  oscillation. 

4 .  Model  Decoupling 

Analysis  of  the  error  equations  provide; 
e,  =*/^<U.V>  -  h,<U, V>  =  U-U 
e  i  =  ><-<U,V)  -  h i  (U ,  V >  =  V(X-3U> -VCX-5U) 

=  X(V-V) -S<UV-UV) 

The  block  diagram  in  Figure  6.4  indicates  the 
possibility  of  decoupling  the  observer  into  two  first  order 
models  by  letting  J^i  equal  zero  <J,i  is  already  zero)  . 

Figure  6.5  shows  a  block  diagram  of  the  discrete 
subsystem.  The  discrete  system  of  Figure  6.6  is  certainly 
stable  if  it's  eigenvalue  lies  within  the  unit  circle. 

|l-a|  <  1*  0  <  a  <  2*0  <  K.  At  <  2 

Under  this  condition  the  system  has  a  steady-state 
error  of  zero  for  a  constant  input  u. 

Hence  the  system  acts  as  positioning  servo-mechanism 
and  is  capable  of  tracking  the  relatively  slow  varying 
inputs  with  very  small  errors. 

5.  Subsystem  Two 

A 

Assuming  that  subsystem  1  is  tracking  u  C  u<t>  = 
u<t>>  the  error  e t  can  be  approximated  as: 

♦  <  A 

Si-  CX  -  SU>  <V  -  V) 

The  observer  is  reduced  to  the  gradient  algorithm 


where  the  gradient  is: 


COftftcSPOMDWVG  tH$,CRET£  MODEL 


^xKt|  =  (  I  -h,  At)--UK  ■+■  At  UH 
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Figure  6.5  Subsystem  1 


Figure  6.7  illustrates  the  observer.  This  result  is 
also  in  accordance  with  the  method  proposed  in  CRef .  6J . 

*  ‘  'V 

Defining  ^ =  (X-SU)V  as  a  new  measurement  the 

*  . 

corresponding  function  h  is  h=<X-5U)V.  The  Jacobian  is 
J  =  J  l:/<Jv=  X-5U. 

In  order  to  simplify  the  stability  analysis  the 
Jacobian,  J,  is  replaced  with  the  signum  function, 
signum(J).  This  is  sufficient  to  guarantee  the  negative 
feedback  of  subsystem  2.  Subsystem  2  is  analogous  to 
subsystem  1  and  will  be  stable  if  0  <  K*.  At  <  2. 
Consequently  this  system  also  performs  as  positioning  servo¬ 
mechanism  and  will  track  relatively  slow  varying  inputs  V(t> 
with  very  small  error. 

The  presented  analysis  leads  to  a  definition  of  the 
decoupled  observer.  This  observer  is  depicted  in  Table  22. 
Appendix  D  contains  the  computer  simulation  results. 

Figure  6.8  is  essentially  the  same  as  Figure  6.4 
with  J*i  =  0  and  Jxi  =  signum<J2e  ).  Therefore  it  is 

possible  to  consider  this  algorithm  a  simplification  of  the 
algorithm  developed  in  CRef.  7].  However,  it  should  be 
pointed  out  that  the  decoupled  observer  can  be  derived 
directly  through  analysis  of  the  two  linear  subsystems 


presented  above 


TABLE  22 


DISCRETIZATION 


-  V<X  -  *U>  -  -  h  z  <U, V) 


V  =  K^SygnCJzi  Je, 


After  discretization  this  becomes 


+  K,ate;K 

V*,,  =  V*  ♦  Kj  At  sygn(Ji2  >e.  K 


Figure  6.7  Subsystem  £ 
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Figure  6. Q  Decoupled  Observer 
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In  the  decoupled  model  it  is  easier  to  control  gains 
since  the  J element  does  not  affect  stability.  It  is 
noted  that  slightly  more  oscillatory  behavior  is  present 
at  the  beginning  of  the  simulation.  This  Is  not  unusual  for 
gradient  type  algorithms. 

Satisfactory  results  for  the  product  K&t  =0.1  are 
expected  if  the  observer  is  analyzed  from  the  sampling 
theorem  point  of  view.  This  presupposes  that  both 
continuous  subsystems  posses  only  one  time  constant  ^  =  1/K. 
The  sampling  theorem  requires  that  t  <  T'/ 2  and  4tK  <  0.5. 

Experience  has  shown  that  the  best  results  for  the 
one-step  integration  are  obtained  if  the  integration  step 
(sampling  time)  is  ten  times  less  that  the  smaller  time 
constant  in  the  system.  This  means  xt  *>  0.1 (1/K)  and 
*tK  *  0.1. 

The  sampling  theorem  is  violated  if  AtK=l.  As  a 
consequence  the  observer  requires  much  more  time  to  extract 
sufficient  information  from  the  measurement  and  to  start 
state  tracking.  The  simulation  results  demonstrate  that 
when  the  product  4tK  equals  unity  convergence  starts  after 
fifty  seconds.  This  is  approximately  ten  times  slower  from 
previous  experiments  where  the  sampling  theorem  is 
satisfied . 

In  the  case  of  4tK  =  0.01  the  system  pole  is  very 
close  to  the  unit  circle.  Therefore,  convergence  is 
relatively  slow. 
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As  demonstrated  by  the  various  computer  generated 
graphs,  the  simulation  results  for  this  system  are 
satisfactory.  Therefore  the  selected  coordinate  system  is 


proven  acceptable 


VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  objective  of  this  thesis  is  twofold.  The  first  part 
of  the  thesis  involves  finding  a  suitable  coordinate  system 
in  which  to  establish  a  satellite  tracking  model.  The 
second  part  o±'  the  thesis  uses  this  coordinate  system  to 
obtain  an  observability  analysis  and  design  an  observer. 

This  latter  portion  is  the  main  thrust  of  the  thesis. 

The  satellite  model  is  based  on  the  coordinate  system 
developed  in  chapter  three.  The  model  13  simplified  yet  not 
unrealistic.  The  main  limiting  assumption  coniines  ooth  the 
satellite  and  target  to  a  constant  speed  ana  heading.  <-t n 
unchanging  direction  is  important  due  to  the  desire  to 
maintain  a  constant  angle  a.  These  restrictions  imposed 
on  the  satellite  are  explained  a  detail  in  chapter  three. 
Satellites  and  targets  (airc  .ft)  would  normally  travel  a 
great  circle  path  which  requires  a  constant  beanna. 
Therefore,  the  assumption  is  valid.  This  basic  premise 
resulted  in  the  spherical  triangle  that  forms  the  basis  of 
the  system  presented  in  Table  5.  Once  the  basic  mocel 
became  established  an  observability  analysis  was 
accomplished  in  chapter  four.  The  basic  observability 
theories  developed  in  CRef.  6]  (for  mechanical  springs)  has 
been  successfully  applied  to  satellite  observability.  The 
system  proved  to  be  observable  in  the  wide  sense. 


Further  studies  in  the  area  of  satellite  tracking  and 
observability  are  highly  recommended.  A  more  realistic 
satellite  model  could  be  devised  in  these  studies.  To 
improve  the  accuracy  of  the  satellite  motion  equations  the 
earth  must  no  longer  be  assumed  spherical  and  the 
satellite's  orbit  must  no  longer  be  confined  to  simple 
circular  orbit.  This  leads  to  a  different  coordinate  system 
by  necessity.  Relative  coordinate  systems  should  be 
explored  in  more  detail.  Appendix  B  discusses  this  brieflv. 
A  possible  approach  involves  using  the  Euler  equations  of 
motion.  A  target  position  with  respect  to  a  satellite  Dnme 
axis  -could  be  transformed  via  pitch,  yaw  and  roll  to  a 
satellite  normal  axis.  Enhanced  satellite  tracking  studies 
should  include  a  target  that  is  not  confined  to  a  constant 
heading.  The  target  should  be  allowed  to  alter  course  to 
further  complicate  the  coordinate  system  and  observer  mode.. 
Both  the  satellite  and  target  should  be  subject  to  random 
disturbances.  This  leads  to  the  development  of  stochastic 
models.  References  6  through  8  are  highly  recommended  to 
anyone  pursuing  this  course  of  research. 

This  thesis  has  successfully  developed  a  model  for 
satellite  tracking.  The  system  has  been  determined  to  be 
observable  and  an  observer  has  been  designed.  The  observer 
model  has  been  successfully  simulated. 


APPENDIX  A 


ELEMENTS  OF  ORBIT 

Satellites  travel  in  elliptical  orbits.  Figure  A.l 
depicts  a  basic  ellipse.  The  distances  a,b  and  c  are 
related  by  equation  A.l. 

c'  =  a1  -  b*  C A.  1 ) 

In  this  equation:  a  =  semi-ma]or  axis 

b  =  semi-minor  axis 
c  =  distance  between  foci 
Eccentricity  is  also  a  parameter  often  used  in 
connection  with  elliptical  orbits.  Eccentricity,  e,  is  a 
meassure  of  the  elongation  of  the  ellipse.  Equations  A. 2 
and  A. 3  relate  eccentricity  to  previously  defined 
parameters . 

e'  =  1  -  (b/a>  (A. 2) 

c  =  ae  ( A  .  3  ) 

Other  elements  include  apogee  and  perogee.  Apogee  is 
that  point  on  the  orbital  ellipse  where  the  satellite  is 
furthest  from  the  earth  (focal  point).  The  perigee  is  that 
point  where  the  satellite  passes  closest  to  the  earth. 

These  new  parameters  are  related  to  a  and  e  by  equations  A. 4 
and  A . 5 . 

apogee  distance  -  d,  =  a(l*e)  (A. 4) 


perogee  distance 


d 


a ( 1 -e ) 


(A. 5) 


Any  two  of  the  above  six  parameters  can  define  an 
ellipse.  For  the  special  case  of  a  circle;  e  =  0,  d  =  d 
c  =  0  and  a  =  b. 

To  describe  satellite  motion  more  information  is  needed. 
The  inclination  angle,  i*,  is  the  angle  of  intersection 
between  the  orbital  plane  and  the  earth's  equatorial  plane. 
Figure  A. 2  illustrates.  Information  on  the  period  of 
revolution  of  the  satellite  is  needed  as  well  as  knowledge 
of  the  precise  time  that  the  satellite  passes  it's  apogee  or 
perogee . 

Some  references  use  information  on  the  ascending  node  to 
determine  satellite  paths  and  positions.  Figures  A. 3  and 
A. 4  illustrate  the  relationships.  [Ref.  5]  uses  the 
parameters  in  table  A.l  as  the  necessary  elements  of  orbit. 

TABLE  A.l 

T  -  period  (in  minutes)  w  -  argument  of  perogee 


i  -  angle  of  inclination  e  -  eccentricity 

A SSP  longitude  at  perogee  t  -  time  at  perogee 

The  elements  listed  in  table  A.l  completely  describe  the 
theory  of  satellite  motion.  Perturbations  are  ignored.  In 
practice  the  motion  of  a  satellite  is  very  nearly  described 
by  these  elements  of  orbit. 
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Figure  A. 3  Elliptical  Orbit  Diagram 


APPENDIX  B 


RELATIVE  MOTION 

The  general  theory  of  relative  motion  is  developed  using 
simple  frames  of  reference.  An  object  in  frame  A  moves  with 
respect  to  that  frame.  Frame  A  moves  with  respect  to  the 
fixed  frame  B.  These  frames  are  considered  to  be  coordinate 
systems  or  sets  of  reference  axes.  Figure  B.l  depicts  an 
object  in  the  prime  system.  The  X-Y  axes  represent  the 
fixed  system.  In  figure  B.2  the  prime  system  moves  with 
respect  to  the  fixed  system  a  distance  d.  The  object  in  the 
prime  system  moves  a  distance  d'  within  its  system. 

Therefore  the  object  moves  a  distance  dTrx.+  L  =  d  *  d'  with  a 
velocity  VrtT*L  =  V  *  V' . 

Further  studies  concerned  with  satellite  tracking  should 
involve  relative  motion.  The  fixed  coordinate  sytem  is  a 
three  dimensional  system  with  its  origin  at  the  center  of 
the  earth.  The  center  of  mass  of  the  satellite  is  the 
origin  for  the  relative  or  prime  coordinate  system.  The 
prime  system  orbits  the  fixed  system. 

The  target  is  detected  by  the  satellite  and  its  motion 
is  therefore  meassured  with  respect  to  the  satellite's  prime 
coordinate  system.  These  relative  coordinates  can  be 


A  .* 


Figure  B. 1  Two  Frames  of  Reference 


APPENDIX  C 


ORIGINAL  MODEL  SIMULATION  RESULTS 

This  appendix  contains  the  simulation  for  the  original 
model.  These  graphs  depict  the  states  X  and  $.  Figure  6.4 
presents  the  block  diagram  of  the  model  simulated.  The 
simulations  are  run  for  10  seconds  and  500  seconds.  This 
insures  that  a  variety  of  suitable  graphs  are  available  for 
analysis.  In  order  to  accompany  test  observer  behavior 
several  different  gains  and  time  increments  are  used  in  the 
simulations.  Specific  combinations  of  gain  and  time 
increment  produce  excellent  results  (as  explained  in 
chapter  6) . 

The  observer  model's  ability  to  track  a  constant  is 
recorded  on  eight  graphs.  The  first  two  graphs  (pages  112- 
113)  depict  state  X  for  constant  angular  velocities  and 
different  gain  values.  The  next  two  graphs  (pages  114-115) 
depict  state  $  for  the  same  gains  mentioned  above.  A  gain 
(k)  of  ten  is  clearly  beneficial  for  both  X  and  $,  since 
then  ,3tk  =  0.1. 

The  following  two  graphs  (pages  116-117)  show  two 
different  time  increments  for  state  X,  constant  gain  and 
constant  angular  velocities.  The  graphs  on  pages  118  and 
119  display  the  same  two  time  increments  for  state  S.  As 


expected  the  smaller  of  the  two  time  increments  (0.01 
seconds)  gives  a  slightly  better  performance. 

The  last  six  graphs  exhibit  the  model's  time  response  to 
sinusoidal  angular  velocities.  These  are  the  results  of  the 
test  used  to  determine  the  model's  ability  to  track  time 
varying  functions.  The  first  three  of  these  graphs  (pages 
120-122)  illustrate  a  target  angular  velocity  of 
•001*cos(Dt)  and  a  satellite  angular  velocity  of 
.0009-cos(Dt) .  The  variable  D  is  varied  in  order  to  test 
observer  response  to  different  rates  of  change  in  angular 
velocities.  The  final  graphs  (pages  123-125)  depict  X  and  S 
with  a  target  angular  velocity  of  .005-cos(Dt)  and  a 
satellite  angular  velocity  of  . 00045 • cos ( Dt ) .  The  variable 
D  is  0.01  in  this  simulation.  This  results  In  sinusoids  of 
smaller  amplitudes  and  longer  periods.  The  results  of 
tracking  the  sine  wave  are  very  satisfactory.  It  is  noted 
that  the  relative  error  for  both  states  is  very  small. 
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DECOUPLED  MODEL  SIMULATION  RESULTS 

These  twenty  graphs  resulted  from  the  decoupled  model 
simulations  (as  described  in  chapter  6).  Figure  6.3 
presents  a  block  diagram  of  the  decoupled  model.  The  first 
five  graphs  (pages  128-132)  illustrate  the  time  response  for 
state  X.  The  time  increment ,  ..t,  is  held  constant  at  0.01 
seconds.  Several  different  gains  were  used  on  each  graph. 
The  state  X  was  followed  closely  by  X.  Different  scales  and 
simulation  times  are  plotted  to  produce  a  variety  of 
suitable  graphs.  A  gain  of  ten  consistently  produces  very 
good  results  in  all  five  graphs. 

The  next  group  of  five  computer  graphs  (pages  133-137) 
exhibit  state  $.  As  above,  the  time  increment,  t,  is  kept 
at  0.01  while  the  gain  (k)  is  varied.  The  time  response  for 
the  lower  gains  in  particular  are  very  satisfactory.  The 
state  $  is  closely  approximated  by  $. 

The  third  set  of  five  graphs  (pages  138-142)  represent 
state  X.  This  group  of  computer  graphs  holds  the  gain 
constant  at  one  and  allows  the  time  increment  to  vary.  As 
expected,  smaller  time  increments  (0.1  and  0.01)  result  in 
superior  tracking  capabilities. 

The  last  set  of  five  graphs  (pages  143-147)  show  $  and 
5.  The  gain  is  kept  at  a  constant  value  equal  to  one  and 
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the  time  increment  is  varied.  As  before,  the  smaller  time 
increments  produce  the  best  results. 

In  all  of  the  above  examples,  the  approximated  state  (X 
or  5)  closely  tracks  the  actual  state  (X  or  5) .  These 
simulation  results  are  quite  satisfactory. 
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